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1 . Introduction 


With the advent of large spacecraft, flexibility has become an in- 
creasingly important factor in the system attitude stability. Early de- 
signs of spacecraft were based on rigid body analysis, according to which 
rotational motion is stable if it takes place about the axis of maximum or 
minimum moment of inertia and unstable if the body rotates about the axis 
of intermediate moment of inertia (see, for example. Ref. 1, Sec. 6.7). 

The erratic behavior of the Explorer I, a satellite stabilized about the 
axis of minimum moment of inertia, prompted a re-examination of the rigid 

n 

body assumption. Indeed, Thomson and Reiter were able to attribute the 
behavior of the Explorer I satellite to energy dissipation resulting from 
the vibration of flexible antennas. This conclusion was corroborated by 
Meirovitch.^ References 2 and 3 used the so-called "energy sink" approach. 
Their main conclusion was that a flexible satellite cannot be stabilized 
about the axis of minimum moment of inertia, leaving as stability criterion 
what has come to be known as the "greatest moment of inertia" requirement. 

For a number of years, no significant additional work on the 
stability of flexible spacecraft was perfomed. Some work on cable-connected 
space stations cannot be really considered pertinent. Some investigation 
that can be regarded as being related to flexible spacecraft is that by 
Hooker and Margulies,^ who model a satellite as "a set of n rigid bodies 
interconnected by dissipative elastic joints," and forming so-called 
" topol ogi cal trees . " 

The first serious attempt to treat rigorously the flexibility ef- 
fects on the attitude stability of flexible satellites can be attributed 
to Meirovitch and Nelson.® Reference 5 investigated a satellite with 
elastic appendages by means of an infinitesimal analysis. It appears that 
Ref. 5 uses modal analysis for the first time in conjunction with the 
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stability of flexible spacecraft. At the same time, Nelson and Meirovitch® 
used the Liapunov direct method to investigate the stability of a rigid 
satellite with elastically connected moving parts. The motion of a satel- 
lite consisting of two rigid bodies connected by an elastic structure was 
investigated by Robe and Kane^ by means of an infinitesimal analysis. Simu- 
lating a spacecraft by a set of rigid masses interconnected by massless 
elastic members, Likins® derived the corresponding equations of motion, and 
indicated that a solution can be obtained by modal analysis. Reference 8, 
however, does not produce an algorithm for the solution of the equations. 

Q 

Thermal effects and solar radiation pressure were found by Etkin and Hughes 
to cause the anomalous behavior of spinning satellites with long flexible 
antennas. The flexibility effects on the attitude motion of spacecraft 
were also investigated by Modi and Berenton^® but the validity of their 
analysis is in doubt, as they restrict the satellite vibration to planar. 

An interesting paper by Newton and Farrell^^ presents a method for 
evaluating the natural frequencies of a flexible gravity-gradient stabi- 
lized satellite. In the process. Reference 11 linearizes the equations of 
motion about the deformed equilibrium. As generalized coordinates, the 

m 

investigators consider complete deformation patterns of the satellite. 

This procedure is not only unorthodox but also tends to limit the number of 

degrees of freedom of the simulation, not to mention the fact that one must 

guess in advance configuration patterns. Moreover, there is some question 

as to the evaluation of the equilibrium configuration. Nevertheless, the 

12 

paper contains some interesting ideas. A paper by Likins and Wirsching 
proposes to introduce the concept of "synthetic modes" in conjunction with 
a "hybrid" coordinate system, where the latter is defined as a set of co- 
ordinates consisting of rotational coordinates of the spacecraft as a whole 
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and modal coordinates for the flexible appendages. This idea, however, 
was introduced earlier in Reference 5. 

All preceding investigations have one thing in common, namely, 
they are all discretization schemes. Some use lumping of the distributed 
mass of the elastic members, a procedure referred to sometimes as spatial 
discretization, and others use series truncation in conjunction with modal 
analysis. In a first attempt to apply Liapunov's direct method to hybrid 
systems from the area of satellite dynamics, i.e., without using any dis- 
cretization scheme, Meirovitch^^’^^’^^ studied the^^stabil ity of spinning 
rigid bodies with elastic appendages. It should be pointed out that the 
term "hybrid" refers here to a system defined by sets of both ordinary and 
partial differential equations, a concept different from that used by 
Likins and Wirsching. Several new ideas were introduced in Ref. 13, such 
as the use of the bounding properties of Rayleigh's quotient to eliminate 
spatial derivatives from the problem formulation and the use of testing 
density functions. 

The ideas of Refs. 13-15 have been pursued by Meirovitch and 
Calico^^’^^ for the case in which testing density functions cannot be de- 
fined readily. References 16 and 17 develop the so-called "method of 
integral coordinates," whereby certain integrals are identified and defined 
as generalized coordinate's. Then, using the bounding properties of 
Rayleigh's quotient as well as certain Schwarz's inequalities for functions, 
a function k bounding the Hamiltonian H from below is obtained, < £ H, so 
that K can be used as a testing function in conjunction with Liapunov's 
direct method. The method of integral coordinates is basically a dis- 
cretization scheme. 

One problem that has received little attention in the technical 
literature is that of deformed equilibrium, which can be referred to 
mathematically as "nontrivial equilibrium." Such problems arise in the 
case of gravity-gradient or spin-stabilized satellites with very flexible 
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appendages that are not aligned with the satellite's principal axes. Find- 
ing the equilibrium configuration can be quite a problem in itself, 

especially if the governing equations are nonlinear. Addressing himself 

18 

to this problem, Flatley obtained the nonlinear equilibrium configuration 
of the Radio Astrononiy Explorer (RAE) satellite by means of an analogue 
computer. Deformed equilibrium has also been considered in Ref. 11, but 
the details are not clear and no plot of the deformed equilibrium is shown. 
In seeking stability statements for the RAE/B satellite, Meirovitch^^ ob- 
tained as a by-product the nonlinear deformed equilibrium, thus confirming 
the results of Ref. 18. 

The present study is concerned with the stability of a hybrid dy- 
namical system about nontrivial equilibrium. It contains many of the for- 
mulations and results of Ref. 19. Qualitative stability statements are ob- 
tained for the RAE/B satellite by both the Liapunov direct method and by an 
infinitesimal analysis. In connection with the infinitesimal analysis, the 
natural frequencies of oscillation about the nonlinear nontrivial equilib- 
rium were obtained by a method developed by the first author of this re- 
20 

port. The method of Ref. 20 considers a state vector consisting of 
generalized coordinates and velocities, where the coordinates include both 
rotations and elastic deformations, and develops an eigenvalue problem in 
terms of real quantities alone. The stability statements of Ref. 19 and 
corresponding statements obtained from the solution of the eigenvalue prob- 
lem agree completely. 

2. Problem Formulation 

We shall be concerned with the motion of a body consisting of n + 1 
parts, of which one part is rigid and n parts are elastic. The domain of 
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extension of the rigid part is denoted by Dg and those of. the elastic parts 
when in undefomied state by (i = l,2,...,n) (see Fig. 1). Correspond- 
ingly, the masses associated with the domains are denoted by m.j 
(i = 0,l,...,n), so that the total mass is m = .E^ m^. The elastic domains 
are rigidly attached to Dq and have common boundaries only with Dg. 

The body m is assumed to move in a central -force gravitational 
field, with its mass center describing a given orbit about the center of 
force C.F.., where the latter is assumed to be fixed in an inertial space. 

In describing the motion of m it will prove convenient to identify 
a system of axes xyz (see Fig. 1) with the undeformed state. The origin c 
of xyz is taken to coincide with the mass center of m in the undeformed 
state and axes xyz themselves coincide with the principal axes of m in the 
same state. Note that the system xyz is embedded in the rigid part Dg but 
is not necessarily a set of principal axes for that part. We shall assume 
here that the nature of the elastic motion is such that the mass center of 
the entire system remains at the origin of xyz. In measuring elastic de- 
formations, we consider reference frames x^y^z^ fixed relative to the 
elastic domains D^- (i = l,2,...,n), where the direction of these axes is 
chosen parallel to that of the elastic deformations. The origin of axes 
x^-y^.z^- is denoted by 0^ and in general it need not coincide with c. 

Next let us denote the radius vector from the mass center c to a 
point in the domain Dj (i = 0 , 1 ,..., n) by i* + r;^ , where the point coincides 
with the position of an element of mass dm^- when the body is in undeformed 
state. The constant-magnitude vector denotes the radius vector from c 
to 0-j ; clearly hg = 0. On the other hand, r^ is the radius vector from 0^- 
to the point in question, and its components represent the independent 
spatial variables associated with a point in the domain D^-. Denoting by i^i*, 
, and the unit vectors along axes x,- , y^ and z^ , respectively, we can 
write h^. + r ^ = (h^^,- + x^Oii * ^^yi ^^zi ^i^lLi ' 0,1,..., n). 
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In describing the elastic deformations, we can use the Lagrangian or the 
Eulerian approach. According to the Lagrangian approach the independent 
variables are those of the body in undeformed state, whereas in the 
Eulerian description of motion the independent variables are those of the 
body in deformed shape. For infinitesimally small deformations the two 
points of view coalesce, but for large deformations they do not. When it 
is necessary to calculate the stresses in a body undergoing large deforma- 
tions, the Eulerian approach is more convenient. Although we shall be con- 
cerned with relatively large deformations, we have no interest in the in- 
ternal stress distribution, and because of kinematical considerations we 
shall find it more convenient to use the Lagrangian approach. Hence, de- 
noting by the elastic displacement vector of dm^- , and recognizing that 
the vector depends both on spatial position and time, we can write it in 
the form u^ = u^ (x^ ,y^ ,z^ ,t)i^ + (x^ ,y^ ,z^ ,t)j^ + (x^ ,y^ ,z^ ,t)k^ . where 

u^ , v^ and w^ are displacement components measured along x^ , y^ and 
respectively. If is the radius vector from the center of force C.F. to 
the mass center c, then the position relative to the inertial space of a 
mass element dm^ in deformed state is given by “ Be lli '*'11 Jdi* 

It should be noted that, by the definition of the mass center, 

(io jm, ■ £• 

In view of the above discussion, the kinetic energy can be written as 



+ Hi) • (ii + ii + Hi)‘l'"i ^1) 

where the first term on the right side of Eq. (1) is recognized as the kinetic 
energy of translation of the mass center c and the second one as the kinetic 
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energy due to motion relative to c. Dots denote derivatives with respect to 
time. Denoting by u the angular velocity of the set of axes xyz, hence also 
of the sets x^y^-z^. (i = 1,2,. ...n), and recalling the expression for the 
time derivative of a vector expressed in terms of rotating coordinates, we 
obtain 


h. 


. + r. + u.. = ul + o)x (hi + r.* + u.-) 

■j *.*1 ♦••■I 4M. ^1 •»! 


( 2 ) 


in which u,! = u^i,- + v,-j,- + wk,. is the velocity of dm.- relative to c due 
to elastic effects alone. Introducing Eq. (2) into (1), we arrive at 

T = jSc + • Jd • JS* “ • / (Hi ♦ n * H1) “ “i*-! 

1=1 ■'m2 


i=l 


m2 


ui * u.' dm 2 

-M I 1 


(3) 


where is the inertia dyadic of the body in deformed state about axes 
xyz. 

Equation (3) is most conveniently expressed in matrix form. The 

• • 

matrix forms of the vectors and a are simply {R^.} and {w}, respectively. 

The inertia dyadic J(j and the term on the right side of Eq. (3) require 

n (i) 

further attention. The inertia dyadic can be written as , 

where (i = 0,1 n) is the inertia dyadic associated with domain 

D^- when the corresponding mass is in deformed shape. The superscript 
i indicates that the dyadic is expressed in terms of the base x^y^z^. 

This would require that we express u in the same base. It is simpler, 
however, to express every in the base xyz instead. Denoting the 
vector r,. by r!^^ and r!^^ when expressed in the base xyz and x.y.z. , 
respectively, and by {rj^h and {r!^^} the associated column matrices. 
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the relation between the two can be written as 

where is the matrix of direction cosines between axes and 

xyz. In a similar fashion, if we denote by and the inertia 
dyadics when expressed in the base xyz and x^y^z^, respectively, and 
t>y the associated inertia matrices, then the relation 

between the two can be shown to have the form 

With the understanding that the inertia matrices imply the body in deformed 
shape, we can drop the subscript d. Moreover, we shall drop the super- 
script i when it agrees with the subscript. Hence, the inertia matrix for 

the entire body, expressed in the base xyz, takes the form = 

n j 

^Iq [^i] We note that [iIq] = [1], where [1] is the unit 

matrix. A similar analysis can be performed with regard to the third 
term on the right side of Eq. (3). It follows that Eq. (3) can be written 
in the matrix form 


1 . T . 1 n T T T 

T = 2 ^^c^ 2 ^^i^ [Ji][ili]{co} + {u} 


[h 


( 0 ) 


mi 


(0) (0) T • 1 f • T • 

+ r. + u. ][^.] {uDdm. + j z {ul} {ui}dm. 
1 1 1 1 1 ^ j 1 1 1 


(4) 


where [h| ^ + r|^^ + u|^^] is a skew-symmetric matrix whose elements satisfy 
.u- u(0) . JO) . (0) 3 ,.(0) . (0) . (0). . ... 


the reletton . u<°) = (h|f t r'°> t u<°>), in which 

^nmJi the epsilon symbol (see Ref. 1, p. 109). Clearly, {ul} represents 
the matrix notation of u^I . 

The potential energy results from two sources, namely, gravity and 
elastic deformations, denoted by Vq and V^, respectively, so that 
V = Vq + Ve- From Ref. 15, we conclude that the gravitational potential 
energy can be written as 
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I^m T TT 

''g " ■ ^ ■ 2r| ([AirCOl][iii] )+^ {Aa> [Ji][Ai]{iia> (5) 

.where tr denotes the trace of a matrix, and {Jig} is the column matrix of 
direction cosines between the direction of the vector^ and axes xyz. 

The elastic potential energy, also known as strain energy, requires 
special attention, particularly in the case of large deformations. No 
general expression, such as for T and Vq, can be written for V^. This is 
so because an explicit form requires the knowledge of the type of elastic 
members involved. For the moment, we shall be content to write 

n 

Vr = E V... (6) 

i=l 

where (i ■ l,2,...,n) is the elastic potential energy associated with 
the member occupying the domain D^* when the member is undeformed. We 
shall return to the elastic potential energy shortly. 

At this point it appears desirable to determine the functional 
dependence of the kinetic and potential energy in order to derive general 
Lagrange's equations of nwtion. To this end, we must specify the nature of 
the elastic members. We shall be concerned with one-dimensional members 
capable of flexure in two orthogonal directions. Any axial displacements 
will be assumed to be a result of change of length caused by the trans- 
verse displacements and not because of axial flexibility. In essence, the 
members are cantilevered bars undergoing large transverse displacements 
(see Fig. 2). Although we shall use nonlinear theory for the elastic 
motion, this will be because geometric nonlinearities and not as a result 
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of nonlinear stress-strain relations . The mass distribution is arbitrary, 
but some of the members carry tip masses. 

Letting the radius vector r,- be aligned with axis x^. when the bar 
is undeformed, we conclude from Fig. 2 that 

t 

.(i) = 


Xi 


r. = X. i . , i = 1 ,2,. . . ,n 
■^1 1 o*! ’ 


and 


u!^^ = u (x. ,t) = V (x. ,t)j. + w.(x. ,t)k. , i = l,2,...,n 

^1 I 1 I II 


(7) 


( 8 ) 


In view of this, the elements of the inertia matrix for the rigid member 
can be written as 


'^011 ^ ^0 ’ *^022 “ ’ *^033 ^ ^0 


*^01 2 " '^021 ^ '^013 * *^031 “ *^023 ' '^032 " ^ 


(9) 


where Aq, Bq, Cq are the principal moments of inertia of the rigid part. 


whereas 

these for member i 

are 








'^ill " J 


("zl 

+ w. 
1 

)^]dx^ + 

miC(hy^ 

+ Vi)2 + 

P*z1 

+ W,)^] 

'^i22 ' j 

f^i 2 

^PiC(hxi+Xi) + 

(h^i 

+ W. 

)^]dx^ + 

miC(hxi 

+ Xi)^ + 

("z1 

+ w,)2] 

*^133 " J 


(hyi 

^ ''i 

)^]dx. + 

mi [(hxi 

+ Xi)^ + 

(hj,, 

+ V,)^] 



r^i 



+ v^.)dx^. 




V,) 


'^il2 '^i21 " " 

Q Pi^^i + 

XlO(byi 

mi (hj^i 

+ x,)(hy, + 

Xi = H 

'^il3 " '^i31 " ■ J 

f)^i 

0 ^ 

XiHhji 

+ w^)dx. 

mi (hj^i 

+ x,)(hj., + 

«1) 

x.=fc 


*i’‘i 


*i“‘i 


x,=t, 
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*^123 '^132 




Pi(hyi + v^)(h 2 i* + w.)dx^. - + v^OChj,,- + w^) 


x.=£,. 


1 = 1 ,2,. . . ,n 


( 10 ) 


Note that and m^. are mass densities and tip masses, respectively, and 
h . , h . , h . denote the coordinates of the points of attachment of the 

A I Jr I I 

booms measured from the mass center along axes (i = l,2,...,n). 

We shall not specify the mass densities at this point. 

The desired equilibrium configuration is that of gravity-gradient 
stabilization. That implies that the mass center c moves in a circular 
orbit with the constant angular velocity (see Fig. 3), and the set of 
axes xyz coincides with a set of orbital axes abc, where ^coincides 
with the direction of the radius vector R^, b is tangent to the orbit 
and in the direction of the motion, and c is normal to the motion. Note 
that the orbital axes abc rotate relative to an inertial space with angular 
velocity ^ about axis c. The orientation of axes xyz with respect to abc 
is given by three angles ej and {oi} depends on these angles and angular 
velocities ej (j = 1,2,3). Because the first term in the kinetic energy, 

Eq. (4), is constant for a circular orbit, it will be ignored in future 
discussions. Moreover, the last term depends on the elastic velocities, 
so that the functional dependence of T is 
• • • 

T ~ T ( 9 j , 6 j , > ^i ’ ^i * '^i ^ » J * 1»2,3 j i — l,2,...,n (11) 

The gravitational potential energy Vg contains the matrix which 

is defined as the matrix of direction cosines between R^. and xyz. Since 
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xyz can be obtained from abc by means of the rotations (j = 1,2,3), it 

v! 

follows that 




J ~ 1,2,3 ^ i ~ l,2,***,n 


(12) 


It remains to establish the functional dependence of V£. This 
requires some elaboration, particularly because of the geometric nonlineari- 
ties involved. First we wish to distinguish between the potential energy 
V£^ due to axial motion, and the potential energy V^g due to flexure. 

Next let us consider Fig. 4 and denote by s^- the distance to any element 
of mass dm^. when measured along the deflected bar and by x^ when measured 
along the original direction of the undeflected bar. We shall assume that 
the bar is inextensional , so that these two distances remain the same, 
s^- = x^ . An element of length along the deflected bar can be obtained 
from 

(ds^.)^ = (dx^ + du,-)^ + (dv^-)^ + (dw^-)^ (13) 


Assuming that du^. is one order of magnitude smaller than dv^. and dw^. 
recalling that ds^- = dx^. , and rearranging Eq. (13), we arrive at 

‘‘“i ' - 2 [(a;;}) * (^) ' = i,2,....n 


(14) 


SO that the axial displacement resulting from the transverse displacements 
is negative. Because for inextensional motion the axial force does 

A I 

not depend on the axial displacement, and, moreover, because a tensile 
force opposes the motion, we have 


TA 


= - z 


i=l 


XI 


j„ ”41 [y * y K 
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where total derivatives have been replaced by partial derivatives in 
recognition of the fact that the displacements are functions not only of 
spatial position but also of time. 

The flexure potential energy is due to the bending displacements v^ 
and w^- . We shall denote the bending moment associated with the displace- 
ment v^- by M^-j and the change in slope corresponding to an element of 
length in deformed state by because they both take place about the 
z^--axis. Accordingly, the analogous quantities associated with flexure 
about y^ are denoted'by and d(()y^ , respectively. It follows that the 
flexure potential energy can be written as 



But the bending moments and can be written in terms of the 
associated flexural stiffness and radii of curvature, as follows 


El 


M, 


yi 




M 


El 


zi 


Zl 


R, 


(17) 


•yi - 'Zl 

where, from Fig. 5, the radii of curvature have the form 


in which 


ds. 


Zl 


R = 


fl + 

'dVj- 

|2- 

L 

[dxj 



Zl 


dxj 


d(j> 


ds 


yi 


ii 


’zi 


Fl X 

fdw,|2i 

= [1 - 

J 


1/2 


dx. 

1 


(18) 


(19) 


Moreover 

dv,- . ‘^'^i - * A 

= tan = tan 


(20) 
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From Eqs. (20), it follows that 


f -1 ^''il 


= d 


r -1 dW,l 

tan T— J-' 

dx^ 


[dx,. 


dx,.^ 


fdvil2 


1 + 


dx^- 


fdw.) 




dx; 


dx^. 


O 

d^w 


1 


1 + 




dx^ 


1 + 


dxi^ 

P dx. 

rdWd^2 1 


"dw j 'I ' 


( 21 ) 


Finally, introducing Eqs. (17) through (21) into Eq. (16), we obtain 

.2 


''eB " 2 1^. ^‘*’zi ^yi ^“^yi ^ ' 2 ^-=1 L^. 


El,, 


(d<l>2^ ) 


zi ds. 


(d^^) 

^^yi 


2n 


1 ^ 


n fij 
i»l ‘o 


EL 



d^v^ 

2 [d2w,' 

CM 


dxi^ 

+ El 


"l+ 

[‘‘XlJ 

>■15/2 ruf*'^i‘ 

J L 

1 1 1 

CJ1 

ro 


Vdx^ 


(22) 


Recalling that the flexural displacements depend also on time, and writing 
binomial expansions for the denominators in Eq. (22), we arrive at 

2 ... . 2, 


n fZ 


'eb 4 J 


'1 


0 


9 t 


El 


Z1 


3X-j2, 


1 


1 f!!i' 

2 ^3X^ 


-r a' 


5 

2 lax^. 


Hx^- 


(23) 
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where the terms involving 9v^/3x^- and 9w^/3x^- are recognized as the cor- 
rections due to the geometric nonlinear effect. 

In view of the above, the potential energy has the general functional 


form 


'^EA ''eB 

= Vp (vl, v", w!, wV) , i = l,2,...,n 


(24) 


where primes indicate differentiations with respect to x^-. 

From Fig. 4, we conclude that we must still account for the dis- 
tributed forces Py.j and p^^-. Regarding these forces as nonconservative, 
and assuming that they do not depend on the elastic deformations, we can 
account for their effect in the form of the nonconservative work 


SO that the total work can be written as 

„ . * M„c " -V * M„c (26) 

where the conservative work has been recognized as being equal to the 
negative of the potential energy. 

The system differential equations of motion, and the appropriate 
boundary conditions, can be obtained from the extended Hamilton's 
principle (see Ref. 1, Sec. 2.7) 

rt« 

(6T + ^dt = 0 (27) 

^^1 
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where all the virtual displacements must be set equal to zero at 

t = t-|,t 2 . Introducing the Lagrangian L * T - V, Eq. (27) can be written 

as 



+ • 0 


(28) 


where the Lagrangian has the functional form 
• • • 

L = L(6j, 6j, v^, v^, v], vV, w^, w^, wl, Wp , j = 1,2,3 ; i = l,2,...,n 

(29) 


It will prove convenient to separate the Lagrangian into that 
associated with the rigid domain Dq and those associated with the elastic 
domains . Hence, let the Lagrangian have the general functional form 
(see Ref. 19) 


L(t) = Lg(t) + E 
^ i=l 


fU. 

L^(x^. ,t)dx^ + L^()l^,t) 


(30) 


where 


Lo(t) = Lo[0j(t), 0j(t)], j = 1,2,3 


L^-(x^-,t) = L^[ej(t),0j(t),v,-(x^ ,t), v,*(x^,t),v-(x^. ,t),vy(x.j,t): 

^i (^•j »^) 3 
i = 1 ,2,. . . ,n 


(31) 


(32) 


L^.(£. ,t) = L^.[0j(t),6j(t),V^(il^ .t),V^(i^,t),W^(t^,t),W^(£^,t)] 


in which Lq is the Lagrangian corresponding to the system in undeformed 
state, Lj the Lagrangian density associated with any point of the elastic 
member i, and L.,- the Lagrangian corresponding to the tip mass. Moreover, 
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represents the length of member i. From Eqs. (30), (31), and (32), we 
conclude that 


3 

6L = E 

j=l 


hi 


dL 


1 " 

n 

+ ^ ■ 

J i=l 

Jo ^ 


" sLi sL-j • sLj , sL-j 

* ITf 357 «“i * ^ 3^^ ‘"1 55f ««1, 


dx. 


8L 3L** 3L 

^ 3Vi-(£i*,t) * 3vY(£:;t) ,t) + 3w,-(^J,t) 


9L, 


3 wV(£:",T)' «Wi-(Ai*,t)j 


(33) 


In addition. 


n 

'Wpc = 

1=1 


0 


(Pyi^Vi + 


(34) 


Inserting Eqs. (33) and (34) into (32), and integrating by parts with 
respect to t, we arrive at Lagrange's equations for the rotational motion 


3L _ 

36j dt 


OL 1 


1.9 6 j 


= 0 , j = 1,2,3 


(35) 


Moreover, integrating by parts with respect to t and x^- , we obtain 
Lagrange's equations for the transverse displacements , and the associated 
boundary conditions, in the form 


A 

ii. 


[3Li] 

_ _a_ 

f3Lil 

*4 


av^* 

at 

[aviJ 

ax^- 

.3vj. 

ax| 



0, 0 < X. < 

1 * l,2,...,n 


and 


(36a) 
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(36b) 


^ - 

a 

pLil 


a 

[9LiV 


axj 

av- 

< ' j 

avj 

at 



/s 

8Li 

5Vi‘ = 0, 6v^- = 0 at x^- 


i = 


■’Li _ 

[’hll 

3vj 3x,* 

av^' 

. • / - 




6v- = 0 at 


0 


Equations similar in structure to Eqs. (36) can be written for w* by 
simply replacing v^. by w^. . 

3. Nontrivial Equilibrium 

Let us consider the case in which Py^* = p^j = 0 and define an 
equilibrium configuration as a set of dependent variables 9j, v^-, w^ 
constant in time and satisfying Lagrange's equations. Because these 
variables do not depend on time, they must satisfy the equations 


and 




avj ax^- 


h’ 


[9Li1 

'J 


avl' 

iJ 


_ 3 

aLi 


av^! ax^- 

avV 

V * J 

avi 


-A. 

aLi 

a 



a^ 

axj 

av^ 



- 0» 0 < x^- < )i^' , i = l,2,...,n 

A 

3Li 

6Vi = 0, 6vl = 0 at x^- = 

i = 1 ,2 , . . . ,n 

aCi 

*^''i = ■ iTv’ 'Sv,* = 0 at X. = 0 


(37) 


(38a) 


(38b) 


as well as a set of equations similar to (38) for w^. We shall denote the 
solutions of Eqs. (37) and (38), together with the set of equations for w. , 
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by 0jQ, v^.q(x^.)j first are constant and the latter 

functions of the spatial variables x^* alone. 

4. Perturbations About Equilibrium. The Variational Equations of Motio n 
The interest lies in the stability of the system in the neighbor- 
hood of the nontrivial solutions Sjq, v^-q(x^-)> shall seek 

stability criteria by means of Liapunov's direct method, and, to this end, 
we let the solutions of Eqs. (35) and (36) and the companion equations to 
(36) have the form 


6j(t) = 9jo + 9ji(t). j = 1.2,3 

v^(x.j,t) = v.j-|(x.|,t), w^(x^,t) = w^q(x^) + w^i(x^,t) 

i = 1 ,2,. . . ,n 


(39) 


where eji(t), v^..|(x^. ,t), w^^(x^,t) are small perturbations. Inserting 
Eqs. (39) into Eq. (30), and expanding a Taylor's series about the non- 
trivial equilibrium, we obtain 


L = L ( 6 jQ , v ^ Q , vjQ , vjQ , w ^ Q , wjQ , w ^' Q ) + 


9L fl + 3L : 


n 

+ z 

i=l 


r 

L^o 


3Li 


3Vi0 


'11 


J. y J. wl + ^^i V" + W + 

3^''n ^ 3^''il ^ 3^''il ^^^1 ^ 


3Vi‘o 


9V!(0 


awio 


* 3W^ *^11 



'3L^. 3L. . 


dx^ + 

J 

3V^ * ••• + ^^*1 

“ ^•i. 


1 3 3 

j=l k=l 


+ 2 .2 i 


3^L 


0.-ie,o + 2 










0 .. 0 . 


sejo’Vo 
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1 


n 




rd 


r^2 


L. 2 *2 2 „2 2 

^ vfi + -^ V?, + 1 vL + —4 V, + w,, 


L^''iO 


2 ’'il ' ^w2 ''il ^ T:;7z ’'ll " gyM2 ’'il ' 3 ^? "il 




^''ib 


lO 


lO 




2 9^L,. 


. •? ^ / d L - 2 

+ — ^ wf, + — -i- w], + -4- wv; + 2 E 


' A. 


“fo ^ ^ •” 


j=l l^®jO^''iO 


— «jl^l 


a^Ci 


32L. 




^ aojoavio + ... + 




f 92 l. 


+ 2 


32u. 


a^L,- 


^av^.oaw.Q ''ilM " av.oa^,.^ " aw.^av^.Q ^l^il " 




A 


v:,vi', + 


.. ..... 'iTii "irii 


^ W!,WV 


5''i0^''i0 


awioawi'o 


dx^. + 


„2 


a^L 


i ‘2 

■ %4^ 


3V 


^ ''il " ''il 


iO 


"''iO 


a^L. 2 .2 3 

+ — 7T- w,.i + — s- W., + 2 I 


f a^L. 


aw 


r il 

iO 


aw 


2 il 
iO 


j=l 


asjoavio 


^®jl''il 


a^L. 


aejoav^-Q 


®jl''il 


• • 


a^L, 

. -r :— 6j,w„ 


aojoawi-o 


+ 2 


a^L^ a^L^. 

avi-oaw^-o ^‘1^1 ^ av.Qa^/.o ''il'^il 




X. = il,. 

1 '■ 


where aL/ae^Q = aL/aej 


(40) 


®j " ®jO’ •• 


, etc. But the term L(8jQ,v.Q,. - . , wVq) 


is constant. Moreover, by virtue of Eqs. (37) and (38) and the companion 
equations for w^. , all the linear terms in the perturbed variables in 
expansion (40) reduce to 


^ al 
j=l ^^jO 


'jl 


n 

+ E 

i=l 


r*^i 


Jo 


at. . 

_1_ V 

av,.„ 


aL. . 


fat, 


dx. + 
1 


av 


iO 


''il 
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"-i • 1 


which are all linear in the generalized velocities ®ji ,w^.-| . In view 
of this, if we retain terms through second order only, the Lagrangian 


becomes 


where 


" "^21 ’’’ll ’’’oi 


1 3 3 g2L . in 

^ j=l k=l 39jO^®kO ^ i=l 


ft. fa^L. .2 9^L.. .2 

I -yL v^, + w. 
0 _3v?g 3*f(, 


3^L4 • • ^ 3^L,- • • d^Li • • 

3V*jQ9W-jQ j=*| ^ ^ ^ 


Vl. p S^L. 2 3 f a^L-. 

r^^-l ^ ^^-l'^il‘^2 ^ T“®jlVil 


3V..rt9W 


iO 10 


j=l [sejQav^.Q 


S^L, . . 1 

* 36„3W„ ®J1“11 

ju 10 ^i ^i 


is quadratic in the generalized velocities. 


^ 1^1 |Jo b=i i 


3 3 9, 


n ri, [3 a'^L. 




. a2i. . a2[. . 

+ i — e.,v., + i— e.,w., + -:; ! — e-iw.i 

aejoav^-O aejoSw^o aej^aw^-O J' 
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J=1 l9e,'o3V 


96j03Vi0 


99j03Wi-o 


80in9W 


3^L. . 92l. . 

+ ^ v.,w., + ^ — v.,w., 

9Vi09Wi-0 ^V^o3W^•o ” ’ * S,^ 


is linear in the generalized velocities, and 


^01 '^1 " 2 . 1 ^ 1^^^ 9ejo30kO ^jl^kl ^i=i Jq 


96jo3W^.O 


9V^03WiO " 


11 2 9^L< 2 ^ f 

X L_w:,w'.', dx. + V., + — i w,, + 2 z — 

jwio3«;'o " "JJ ' L>''?o »“io 


j=l ^30jo9V^O 




Xi = 


is free of generalized velocities. 

In view of the above, the perturbed Lagrangian can be written in 
the general functional term 


L - L ( 0 j 1 ,0 j 1 , v^ 1 , V j 1 , vj 1 , . . . ,wj .| ,w^ 1 ) , j - l,2,3j i - 1,2,. ..,n 



Consequently, the variational equations can be written in the form of 
the Lagrange equations, Eqs. (35) and (36), but with the subscripts j and 
i replaced by jl and il, respectively. Unlike Eqs. (35) and (36), the 
variational equations possess trivial equilibrium. 

5. Discretization by a Rayleigh-Ritz Approach 

The variational equations discussed in the preceding section con- 
stitute a set of hybrid differential equations, in the sense that the 
equations for the rotational motion are ordinary differential equations 
and those for the elastic displacements are partial differential equations, 
where the latter are subject to given boundary conditions. It will prove 
convenient to transform the system into one consisting of ordinary 
differential equations alone. This can be done by using a discretization 
procedure based on the Rayleigh-Ritz approach. Indeed, let us introduce 
the notation 


6j^(t) = qj(t), j = 1,2,3 


Vtt(x^ ,t) = 


p+3 

2 i>i(xi)q.-(t), 
j=4 ^ 


2p+3 

Wii(Xi,t) = il»j(xT)qj(t) 


3p+3 

'^21^^?*^^ ~ ^ ^■j(x2)Q4(t), 

^ j=2p+4 ^ J 


4p+3 

W9, (Xp,t) = I i|«.(x2)q.(t) 

j=3p+4 ^ ^ J 


(47) 


(2n-l )p+3 2np+3 

Vni(Xn»t) = E <l'i(x )q.(t), w ,(x ,t) = E i/ii(x )q.(t) 
j=2(n-l)p+4 J J " j=(2n-l)p+4 ^ 

where <l>j(x^.) and 'i'j(x^*) are admissible functions, taken as the eigen- 
functions of the linearized system. With this notation, Eq. (43) can 
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be written in the matrix form 


T 2 , = 1 {q(t))V]{q(t)) 


(48) 


where [m] is a constant symmetric matrix having the elements 
m^,. = , j,k = 1,2.3 




(49a) 


m 


rti a^L. 

'0^ = /o 


r 






Xj » X, 


j = 1,2,3; k = 2(i-l)p+4, 2(i-l )p+5,. . . ,(2i-l )p+3, 

i = 1 ,2,. . . ,n 


(49b) 


fii a^^L. 
m,-|. = ; — (|). (x,.)dx-. + 


r ^2 


0 30jo3Wio 




r-^*k<xi> 


*®jo’"io 


X^ = £■{ 


j = 1,2,3; k = (2i-l)p+4, (2i-l )p+5,. . . ,2ip+3, i = l,2,...,n 


a^L- a^L-. 


0 


32l.. 


L=''10 


(49c) 


Jx,. = 


j,k = 2(i-l)p+4, 2(i-l)p+5 (2i-l)p+3, i = 1,2 n (49d) 

.Hi ^2| 


In' 4j(xi)*k(x,)<lx, * 4j(x,)4k<xi) 


m 


jk Jq av^o^Wjo 


3Vio9Wi*o 


Xi = t. 


j = 2(i-l)p+4, 2(1-l)p+5 (21-l)p+3. 


k = (21-1 )p+4, (21-l)p+5....,21p+3, 


1 - l,2....,n (49e) 
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j,k = (2i-l)p+4, (2i-l )p+5,. . . ,2ip+3, i = l,2,...,n 


(49f) 


On the other hand, Eq. (44) leads to the matrix form 


= (q(t)}^[f]{q(t)} 


where [f] is a constant square matrix with the elements 


^jk ” 30 . “ 1*2,3 

JK 36jo^®kO 


^ ' 36 ai * L ' e ' ai ' 

0 L^®jO^''iO Jx. 


(51a) 


j = 1,2,3; k = 2(i-l)p+4, 2(i-l)p+5 (2i-l)p+3, i = l,2,...,n 

(51b) 


f^i ^^*-i r ^^^i 

~ „ ] — K(x,-)dx- + — ^ 

Jo 90jo8Wio 36j03Wio 


Xi = s,j 


j = 1,2,3; k = (2i-l)p+4, (2i-l)p+5 2ip+3, i = l,2,...,n 


f -k = j ^ 

JQ 30kO3ViO ^ 


r 32 l. ] 

(x-)dX. + «|>..(X.) 

' [aeko3vio J Jx^. = 


j = 2(i-l)p+4, 2(i-l)p+5,...,(21-l)p+3, i = 1,2 n; 

k = 1,2,3 


(51c) 


(51d) 
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f,-i, = — » '('•(x. )dx. + — ; ’l',-(x,-) 

Jk Jo 36k0^'^i0 ^ L^®kO^'^iO ^ Jy- = a- 


Xi ■ 


'0 ‘’'^kO^’^iO 

j = (2i-l)p+4, (2i-1 )p+5,. . . ,2ip+3, i = l,2,...,n; k = 1 


'jk 


fili 9"k- 


] — <l'j(Xi-)'l'k(xi)dXi + 


9^L. 




0 9V.q3W.q 

j = 2{i-l)p+4, 2(i-l)p+5 (2i-l)r+3, 

k= (2i-l)p+4, (2i-l)p+5,. . . ,2ip+3, 






i = 1 ,2,. . . ,n 


fi-i 9'^L.. 


^jk 


r 92 l.. 


^ ? — <^k(xi )'«',• (xJdx. + — <{>|^(x.)t|».(x-) 

0 9v,.o3w,.o " 1 J 1 1 .3Vio3Wio ^ ^ 3 \ 


j = (2i-l)p+4, (2i-l)p+5,...,2ip+3, 

k = 2(i-l)p+4, 2(i-l)p+5,...,(2i-l)p+3, 
Finally, from Eq. (45), we can write 

■’’Ol - Vi = - ^q(t)>^[k]{q(t)> 


i = 1 ,2,. . . ,n 


where [k] is a constant symmetric matrix with the elements 

A_ 

^jO^®kO 


•^jk " ■ 30,.f,30Ln ’ 


'jk 


= -f 


£i 9^L^. 




92l 


L^V''iO 


b~~ ^k^^i^ 


JXi = £^- 


j = 1,2,3; k = 2(i-l)p+4, 2(i-l )p+5,. . . ,(2i-l )p+3, 

1 l,2,«**,n 


,2.3 (51e) 


(51f) 


(51g) 


(52) 


(53a) 


(53b) 
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~1 

S‘k ~ ~ Jo 36 jo 9 Wio 


■ 92|_. 

90jO3Wio 


Xj " K^-j 


j = 1,2,3; k = (2i-l)p+4, (2i-l )p+5,. . . ,2ip+3, i = l,2,...,n 

(53c) 


*^jk 


0 


92l. 


9^L-j 9^L,- a^L,- 

avfo J ^ ^ 3 v :2 J 1 k 1 3y„2 J 1 K 1 




) 


+ — ; — J- U-(x.)(frj;;(x.) + <!>'.' (x.)<t'Mxi) 

9 v!o 3 V^'o I j 1 ^ J ’ k 1 J 


dx^. 


92l 




2 


X, = Jj 


j,k = 2(i-l)p+4, 2(i-l)p+5,...,(2i-l)p+3, i = l,2,...,n (53d) 


*^jk ■ " 


;0 - 


■ 92l. 


X. = Ai 


j = 2(i-l)p+4, 2(i-l)p+5,...,(2i-l)p+3, 


k = (2i-l)p+4, (2i-l )p+5,... ,2ip+3 


i = 1 ,2,. . . ,n (53e) 


•^jk ■ ■ 


f^i 

0 




3Wio 


iO 


iO 


A. ) 

iW'" 3 w" + <l'j{XiHk(Xi)]K 

iO iO ^ -'J 
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Vl, 


L’"io 


2 ^ 


Xi = 


j,k = (2i-l)p+4, (2i-1 )p+5,. . . ,2ip+3, i = 1,2 n 


:53f) 


Introducing Eqs. (48), (50), and (52) into Eq. (42), we can write 
the Lagrangian in the matrix form 


L = ^q}^[m]{q} + {q)^[f]{q} - ^q}'^[k]{q} 

Using the approach of Ref. 21 (see Sec. 3-4), we can write Lagrange's 
equations in the matrix form 


(5A) 


d |lLl _ /|]^| = ^0} 

dt [aqj [3q'^ 


(55) 


Hence inserting Eq. (54) into (55), we obtain the equations of motion 

[m]{q} +([ff - [f]){q) + [k]{q} = {0} (56) 


so that, introducing the notation 


[g] = [f] - [f] 


(57) 


where [g] is a skew-symmetric matrix, [g] = -[g], we obtain 


[m]{q) + [g]{q} + [k]{q} = {0} 


(58) 


where [m] is identified as the inertia matrix, [g] is a "gyroscopic" 
matrix and [k] is a stiffness matrix which includes terms due to elastic, 
gravitational, and centrifugal effects. 
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6. Liapunov Stability Analysis 

We shall seek criteria for the stability of motion in the neighbor- 
hood of the nontrivial equilibrium by means of the Liapunov direct method. 
This is equivalent to the problem of stability of the perturbed motion 
about the trivial solution. In terms of the discretized system, the 
perturbed motion is described by the vector {q(t)}, so that the interest 
lies in a stability analysis about the trivial equilibrium (q) = (0). 

It was shown in Ref. 15 that the Hamiltonian is a suitable 
Liapunov function for the type of problem at hand. Assuming that the 
system possesses a certain amount of internal damping, however small, the 
equilibrium is asymptotically stable if the Hamiltonian is positive 
definite. In terms of the perturbed variables, the Hamiltonian has the 
form 

^ " ^21 ■ ’’’Ol ''l = ^{q(t)}"^[k]{q(t)} (59) 

But the function T 21 is positive definite in the generalized velocities 
qj(t) by definition. Hence, if the function 

*K = |-{q(t)}‘^[k]{q(t)} (60) 

is positive definite in the generalized coordinates Pj(t), then the 
Hamiltonian is a positive function in the generalized coordinates and 
velocities and the equilibrium is asymptotically stable. The function 
< is positive definite if the matrix [k] is positive definite. Whether 
[k] is positive definite or not can be ascertained by means of Sylvester's 
criterion (Ref. 1, Sec. 6.7). The matrix [k] will be referred to as a 
Hessian matrix. 
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7. Natural Fre uencies of the Complete Structure 

The Liapunov direct method provides qualitative information con- 
cerning the stability or lack of stability of an equilibrium configuration. 
Similar information can be extracted from the system of equations (58) 
via the eigenvalues. In addition, the eigenvalue problem yields results 
of a more quantitative nature in the form of the system natural frequen- 
cies and the normal modes for the complete structure, where the latter are 
defined later. It turns out that Eqs. (58) lead to an eigenvalue problem 
of a special nature. The nature of the eigenvalue problem can be con- 
veniently discussed by converting the set of equations from second order 
to first order. Indeed, if the configuration vector {q(t)> is of 
dimension N, then we can introduce the 2N-dimensional state vector 
{x(t)> in the form 


( 


{x(t)> = j 

I 


{q(t)> 


(61) 


No confusion should arise from denoting the state vector by {x(t)>, because 
the symbol x.j used to denote the position of a point in the elastic 
members represents a spatial coordinate independent of time and not a 
time-dependent generalized coordinate. Accordingly, if we introduce the 
2N X 2N matrices 


[M] 


[im] [or 

[[0] Ck]J ’ 


[ 6 ] 


■ [g] [kr 

.-[k] [or 


( 62 ) 


then the set of N equations (58) can be transformed into a set of 2N 
first-order equations having the matrix form 
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[M]{i(t)} + [G]{x(t)} = {0} 


(63) 


where [M] is symmetric and [G] is skew-symmetric, 

[M] = [M]T , [G] - -[G]T (64) 

because [m] and [k] are symmetric and [g] is skew-symmetric. 

The matrix equation (63) is of the special form treated in Ref. 20, 
so that the eigenvalue problem can be solved by the method developed 
there. Hence, letting 

{x(t)> = e^^{x} (65) 

where X and {x} are constant, we obtain the eigenvalue problem 

x[M]{x} + [G]{x} = {0} (66) 

It is shown in Ref. 20 that the solution of the eigenvalue problem (66) 
consists of 2N eigenvalues x^ and eigenvectors {xl^. (r = 1,2,...,2N), 
where the eigenvalues consist of pairs of pure imaginary complex conjugates, 
Xp = ± io)p, and the eigenvectors also consist of pairs of associated com- 
plex conjugates {xl^ and {x*}^ (r = 1,2,...,N). Moreover, the eigenvectors 
are orthogonal in a certain sense. Reference 20 provides an algorithm where- 
by the eigenvalue problem can be solved in terms of real quantities. The 
method will be used later in this work to solve the eigenvalue problem for 
a specific spacecraft. 

8. Lagrange's Equations in Explicit Form 

Lagrange's equations, Eqs. (35) and (36), are written in a 
general form. Before obtaining the nontrivial equilibrium and the 
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corresponding variational equations, we must express them in a form in 
which the various coordinates appear explicitly. By virtue of the assump- 
tion that the satellite mass center moves in a circular orbit with 
orbital velocity n, we can replace K/R^. by ?r in Eq. (5). Moreover, the 
first terms in Eqs. (4) and (5) can be ignored because they are constant. 
In view of this, if we recall that the Lagrangian can be written as 
L = X - Vq - - V^g, then we can substitute Eqs. (4), (5), (15), and 

(23) into L, and obtain 

L(t) = i + Tj + i a^trCjC’)] 


f - 1 .1^ * w;2) 


* Elzi-vfo - |vl^) t EIj,,wf(l - |w|2)]dx, 


(67) 


where is the axial force at any point of the slender rod, and 


[J«»] = ? = S 

i=0 i=0 


(68a) 






Xi = 


(68b) 


1 " 
' ? 1 ^, 


I' 


*’*'1 • T • • T • 

p^-{u^.} {u.ldx^. + m^.{u^} {u.} 


X. = 


(68c) 
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in which is the inertia matrix of the body in deformed state in 

terms of the reference system xyz, {K} is an angular momentum matrix due 
to the elastic velocities, and is the kinetic energy due to the elastic 
velocities. The elements of [J^.] are given by Eqs. (9) and (10). Intro- 
ducing the notation 


Jill = OiC(hy,-+v^-)^+(h2i-+w^-)^]. Jin(^) = 




Ji22 = PiC(^1xi■^^i)'^■^^^zi■^'^i^^^’ '^i22^^i^ "" m^*[(hxi"*-x^.)^+(hzi+Wi)^] 


Xi-«-i 


'^i33 = PiC(^i+XT)Mhy^-+v^)^]. J^aaC^i) = m,-[(hx,-+Xi) 2 +(hy,.+v.) 2 ] 


X. = £. 
1 ^ 


'^il2 “ *^121 ~ "Pi (^xi'^’^i ) ^^yi''''^i ^ ’ *^il2^*'i^ ~ '^i21^*'i^ ~ “''^i^'^xi 


'^il3 ~ *^i31 “ "Pi (lixi'’'^i^^*^zi‘'^i ^ » *^il3^^i^ " ’^iSl^^-i) " -nt.j(h 


+x^.)(h^^+w.j) 


XI 


X* = il • 


•^123 ° '^132 “ -Pl0'yi+''i)(''z1-«'l)' Jl23<*i> ' ’^132(‘'l) * 


+V^.)(h2i-+w^-) 


i - 1 7 
1 "■ I yC j • • • jfl 


(69) 


we can write 



In a similar way, from the second of Eqs. (68), we have 


n 

{K} = E 
i=l 


/•£ 
■ 0 




(71) 


In view of the above, the Lagrangian densities can be written as follows 

Mxi,t) = i + i j fi^tr[j{°b 


I - I Pxi(v]2 + wl2) . 1 Ei^. . I vi2) 


- ^Ely,. w;'^(l - |wj2), i = l,2,...,n (72) 

whereas the parts of the Lagrangian associated with the discrete masses 
are 

Li(ti.t) = i {o)}T[j{°^(jl^.)]{a)> + {o)}^{K.(t^.)} + 7m.{u.}^{u.} 


- f (73) 

From the context it should be obvious when brackets and braces denote 
matrices in Eqs. (67) through (73) and when they do not. 

Substituting Eq. (67) into Lagrange's equations for the rota- 
tional motion, Eqs. (55), we obtain 


(U)} + {K} - 3 p2 }T| [j(0)3^„ ^ 


isoj -a- 


It [tI^ ~ ^ " ^’2*3 


(74) 
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Moreover, Lagrange's equations for the transverse displacements v^- are 


? <“> * f“>^ (sir + 1- tr [3|'»] 


m 1 j 


r~A0)- 


- si 


4 - {KJ + Pj V,' 
av^- ^ 


_ _a_ 

3X^ 


(- PxiVi +|eI 2,- vi.v-’^) + -^ v! (1 


c f 2 

“ 2 " ^i ) 3 ^yi * 0<Xj<i^ , i“l,2,...,n 


(75a) 


which are subject to the boundary conditions 


■Cxi - I E'zi ) *1 * jl^ [EIz1 v; (I - f Vi2)] 


? 5I7 Cj'"’]'”’ * f ^ tj"”] - 3 <‘a>^ 5I7 to"»](s,n 


- i? sir * »i = 


EI,i vj 0 = 0 


at x^.=£^. , i 


= 1 ,2,. . . ,n 


(75b) 


f ^ — 0, v^ “ 0 atx^ “ 0, 1 ~ l,2,*««,n 


(75c) 


Similarly, Lagrange's equations for the displacements w^. are 


1 r ,T 

j {o)> 


-i- 

aw,- '■ 1 -■ 


{oj} + {(»)} 


aw 


7 ^*^i^ 


+ 1 tr 


_L. 

^awi 1 
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9W 


r 


f'a) 


_9 

9t 


{(I)} {K} + p-jW,- 

9*i ^ ^ 


5l7(-^iwi*fEv,w;w;'V^c-Ev,w, 0 


- f ■'■^21^^’ ^ ^ ^ ^ l, 2 ,...,n 


(76a) 


which are subject to the boundary conditions 


(fx1 - I E'z1 w?^) «i ♦ ^ [E'yi "i't’ ■ 


5 „.2> 


+ in2 


(«)T ^ [jj'»] („) * tr ^ [j(»)] -3 {ta>^ ^ tJ^^bcta) 


2 - 

3t 


4- {K} + Pi-w.. 

9W,- ' ' 


= 0 


Ely, wf(l - I wj^) " 0 


y at X =ii^. , i = 1 ,2,. . . ,n 


(76b) 


w^. = 0, wj = 0 at X. = 0, i = l,2,...,n (76c) 

9. Equilibrium Equations in Explicit Form 

For a gravity-gradient stabilized satellite, the angles Oj(j = 
1,2,3) are measured relative to an orbiting system of axes. The orbit 
being circular, with the orbital angular velocity being equal to n, the 
orbital axes rotate relative to an inertial space with angular velocity 
about an axis normal to the orbital plane. This axis is denoted by c (see 
complete definition later). Hence, the angular velocity matrix {u} can 
be written as 
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{O)} = + {co}^ 


( 77 ) 


where = {^c(0j)} is the matrix of direction cosines between axis c 
and the reference system xyz, and {a)}p = {oo(9j,9j)}^ is a matrix whose 
elements are the angular velocity components of system xyz relative to 
the orbital axes. They are linear combinations of the velocities Oj 
(j = 1,2,3). 

The equilibrium equations can be obtained by deleting from Eqs. 
(74) - (76) all the terms involving derivatives with respect to time. 

This implies that we can replace {u} by in these equations. Hence, 

the nontrivial equilibrium must satisfy the general equations for the 
rotational motion 

(Ac> - ^ Ua> = 0, j = 1,2,3 (78) 

J J 

as well as the boundary- value problems defined by the differential equa- 
tions 


1.2 


^ * tr ^ [j!'”] - ^ 


* afr [Cxi - I EIzi I’i] - ^ C“zi*lO - 1 = 0’ 


0 ^ x^ ^ , 1 ~ 1 ,2 , • . • ,n 


(79a) 


and the boundary conditions 


•Cxi - I Elzivf )v: + [EI,iVf(l - I Vj^)] 1 1 .2 
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,«»T ,T 9 r,(0)- 




- ElziVi'C - I 'if) = 0 


= 0 


(79b) 


v^- = 0, vj = 0 at x^. = 0, i = 1 , 2 ,. . . ,n (79c) 

Moreover, it must satisfy a set of equations similar in structure to Eqs. 
(79), but with Vj replaced by w^. . 

10. The Variational Equations for the Discretized System 

The variational equations for the discretized system were obtained 
earlier in the form (58), where the matrices [m], [g], and [k] are defined 
by Eqs. (49), (51), (53), and (57). Although the equations just mentioned 
have the advantage of revealing the symmetry of [m] and [k] and the skew- 
symmetry of [g], the formulas for deriving the elements of the matrices 
are not the most suitable from a computational point of view. Indeed, we 
wish to present a procedure whereby the actual derivation of the variational 
equations is performed by a digital computer. 

Consistent with earlier notation, we shall denote quantities 
associated with equilibrium by the subscript 0 and perturbed quantities 
by the subscript 1. With this in mind, we can write the Lagrangian in 
the form 


L. = 


L = Lq + 




'eo 


(80) 


(81) 
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and 


+ {o3}|[J^°hi{a)}o + {to}}{K}Q + {o)}|{K}^ + {w}J{K}i + 

+ tr[J^°hi - 3fi^{£a}|[J^°^]Q{i3}0 - f 

- I n2{£^}J[j(0)]^{,^}^ . 3fi2{A3}|[j(0)]^{A3}o - Vet (82) 


in which 


{oj} = 


{0)}q + {u}^ 


3 

0{£c}q 

1=1 


'®i0 


[O)} 


Pi + 


,^®i0 


{O)} 


+ 1 I I 
2 i=l j=l 


LI 
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i09®j0 


■{oj} 


q,qj + 2 




■{cj} 


‘^j^i 


(83) 




30^0 


Oa) 


1 3 3 

+ r E I 


f 


2 lil 


qi<ii 


(84) 


[J<°>3 = " [tifCJJUi] = " 

i=0 ^ ^ i=0 


'f‘'[j,]dx, + [Ji(4,)] [t,] 

0 j 


= ♦ I [‘1^^ 


1=0 JO 



»Ji 1 
3Vi0 


'11 
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3 J -j 

^ 1 

” 2 ' 
a^J^- 

2 

a^J^ 

»» II X ^ 

a^J^ 

2 

Wii 

9WiO 



Vil ^ 

) 

_av^.o9Wio_ 

Vi^Wii + 2 

O 

<r^ 

1 


dx, 


9 ^ 1 ( ) 


av 


iO 




9 ^ 1 { ^ T ) 


1 


L 


aw 


iO 


Wii (J.-i ) +2 


1 ia^J.j(i,^) 


av ^ 
iO 


v,^M) 




avio9w^-o 


v ^1 + 2 


a^J^.(i^.) 
7 


3^0 


w. 2 (^i) 


[Ml 


(85) 


n 

{K} = z 
i=l 


r 


{K^-}dx. + {Ki(ili)} 


{K)o,= 0 


{K}^ = 


n 

z 

1=1 Jo 


f 9 K^ 1 
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iO 


'i 1 " 


f 9 K^ 1 


^ 

awio 


faK.. 


'il 


f 9 K, 
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iO 
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’“iol 
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av^-oSw^.Q 








y vj^w^-^ 


dx^- + i 
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iO 
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aK.()i.) 

-- 1 1 
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iO 










9 K^(a^) 


awi 


iO 




3 ^K.(Jl.) 




Vji («,-j )w|i (£j ) 


+ 


^ 9 ^K.(£.) 


av]o9w^.o 


vji(Ai)w^l(«.^) (86) 
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Note that Eqs. (83) - (86) represent Taylor's series expansions of 
("a>» and {K} about equilibrium. Moreover, we have 

Teo = 0 (87) 


and 
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(90) 


To obtain the variational equations in terms of the discrete co- 
ordinates qj(t) (j = 1 ,2,. . . ,2np+3) , we must insert the modal expansions 
(47) into L-| and perform the indicated integrations over the spatial 
variables x^- (i = 1,2,. ..,n). Because the resulting expressions are very 
lengthy, we shall not write them explicity, but proceed with the derivation 
of Lagrange's equations instead. To this end, it will prove convenient to 
denote constant terms by the subscript q and terms that are linear in the 
generalized coordinates q,-(t) and generalized velocities q^(t) by the sub- 
script i. This enables us to write 


3qj 








3 T 
^ (u)}^ 


3qj 






j = 1.2.3 


(91a) 
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j “ 4,5,. . ,2np+3 


which enables us to write Lagrange's equations for the perturbed motion 
in the compact form 
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8 r 1 T 
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3q, 
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(O)}, 
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'El 


j = 4,5,. . . ,2np+3 


(92a) 


0 

(92b) 
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As mentioned already, the advantage of Lagrange's equations (92) over 
those derived in Sec. 4 is that Eqs. (92) permit automatic derivation by 
means of a digital computer. 

Before specializing the equations to a particular satellite, let 
us derive an expression for the axial force P^.j in terms of matrix nota- 
tion. The axial force is due to centrifugal and differential gravity 
effects. Introducing the modified potential energy density associated 
with member i 

V, mod = - + tr[3|“>] - 

- 3{£3>'''[j5'”(xp]{ii3))s(xi - (93) 


where the terms inside parentheses and multiplying 6(x^.-t^-) are due to the 
tip masses, the axial force density can be written in the form 


Pxi'Xi) 


^'^i mod 


9X, 


" )T[3!°>]'{£^) + trCj!*”]' - 3{43}^[3!'”]'{ta)) (94) 


in which we introduced the notation 


[3<»>]' =^[3!">]*^[3S°>(x,)]«(x,-t,) 


(95) 


* 2 * 2 

Observing from Eq. (67) that P^.,* is multiplied by (v.| + w^l ), we ignore 
transverse terms in so that using the first of Eqs. (68) and 


any 


Eqs. (69) we obtain the approximation 
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[j''”]' » 2(h^i + Xj) /p, * [t,f 


0 0 0 
0 1 0 
0 0 1 


Ui3 


Inserting Eqs. (95) and (96) into Eq. (94), we can write the axial 
at any point x^- in the form of the integral 


" 9 T f^i (0)' 

Pxi = ^ t Ji (Ci)dc,.]{iic} 


+ tr[ 


rii - 


jj°>'(?^)d5,] - 3{«,)^[ P jj“’'({,)d5, •]{«,}) 

X • *' Xj 


Tr ;(0)' 

Xd 


where, assuming that p^- = const, we have 



(5j)d5^] = ^PiC(hxi+^i)^ 


Olx,+X,)^] 



0 Oj 

1 0| 

0 1 ’ 


It follows that the desired expression has the form 
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force 


(97) 
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0 
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( 99 ) 
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0 

1 
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0 

1. 



n . The RAE/B Satellite. General Formulation , 
a. Equations of motion 

Next let us specialize the equations to the case of a satellite 
consisting of a rigid core with six flexible booms, as shown in Fig. 6. 
First, we wish to determine the matrices [£,i] of the direction cosines 
between axes and xyz. From Fig. 6, it is easy to verify that 
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S6_ 


where sa = sin ct, ca = cos a, S6 = sin 6, and c6 = cos s. Moreover, to 
write the angular velocity matrix {u} in explicit form, we must specify 
the rotations e- (j = 1 , 2 , 3 ). Assuming that system xyz is obtained from 

J 

system abc by means of the rotations 02 about y, -0i about x, and 03 
about z, and recalling that axes abc rotate about c with the constant 
angular velocity i::, matrix {w} can be shown to have the expression 
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ce-| S 03 Ol 
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^1 

S 02 $03 - SO] C 02 C 03 ^ 
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C9-] C 03 0 -i 
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L 0 

S0-, ij 

1 0 o 1 

r 

1 = sin =1 , ce-| = cos 0 ] , 

etc 

Because the direction 


( 101 ) 


radius vector coincides with that of axis a at all times, the direction 
matrix {£3} can be written as 


{£,} = 


C02 

-(c02 S 03 + S9i S02 C 03 ) 
C0-| S02 


(102) 


It will prove convenient to rewrite matrices {03} and {£g} as follows 

M = [933 [e*]-, {0} + {0)3 + n [9]3[0]-|{0}2 003) 

where 
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(104) 


Introducing Eq. (103) into- (4), and recalling Eq. (68), the kinetic energy 
becomes 
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T = ^ {-4 + {oa}''’ {K} + Tf 

= i [°*]j [°]J [e]3 [9*], {M + 0 )}^ [o*]| [9]T [j(0)] (0)3 

+ 1 i^}] -"'3 + [e*]i [6]J {K} + (K} + T^ 

+ ■ [v]| [^l3 [e]3 [9*]^ {«} + ■: ^ [-]{ [9]T [j(0)] {.-,}3 

^ [“]{ [e]^ {K} + fejT [-q^T [-j(O)^ 

= 1 { 9 }T [ 9 *]| [J*] [e*]] {§} + {e}V]j [e]J [J] { 9)3 

+ {9}T [e*]| {K*} + 1 J 33 + K 3 ^ + 

+ {e}J [e]| [J*] [ 9 *]^ {0} + {6)2 [e]j [e]^ [J] {0)3 

+ n {0}J [e]} {K*} + 1 {0}T [e]T [-J*] j-g]^ ^g}^ (105) 

where [J*] = [e]^ [63] and {K*} = [e]^ {K}. Moreover, inserting 

Eq. ( 103 ) into ( 5 ), and recognizing that {^13}”^ = - [q]]" (0^3* 

we obtain the gravitational potential energy in the form 

V5 . - ^ a2 tr + I ■(-(e)^'^ [e]{ [e]J [jf"): [ejj [e], (-{e^ 

= - in^ tr [o(°)] * 3 a 2 ({e)^)’' [ 9 ]j [J*] [e], {e}^ ( 106 ) 

where primes indicate differentiation with respect to 02- Expanding the 
matrix involved in T and Vq, and recalling that L = T - Vg - - V^g, the 

Lagrangian L can be written in the form 

L = i [J*.| ef + (J22 c^e-| + J33 s^e-| + 2 J23 seicei) 
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(107) 


To obtain Lagrange's equations for the rotational motion, we intro- 
duoe Eg. ( 107 ) into Eqs. ( 35 ), and obtain 


[(J33 - J22) S01C01 + J23 (c^ei - s^0i)]02 + (J*2 SO1 - J*3 C0i)0^02 
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+ (-J22 so*! + 00*1)0262 ~ (K2 50*1 - Kj c0*j)02 

+ Q ^^^12 *^13 ^ ^'^12 ” '^IS 

■*■ ^*^33 ■ " s^ei)c02 - 4 J 23 ceiseic62]92 

^ if * ^c 

- (0^3 S61C02 + J23 C0ic02)e3 - K2 C0-|C02 - K3 S0-| C02> 

+ [4 J*2 C0'|S02C02 "•■ (^22 " '^33^ S0]C0] (c^02 ■ 3s^02) 

+ 4 J*3 S0]S02C02 + O23 “ s^ 0 i)( 3 s ^02 - c^02)3 

“ ^ [Jn 0] ” (Ji2 ^ *^13 ■ *^13 ®3 “ 

+ n (J-|i S02 + J-|2 S01C02 - 0-I3 ce-)C02)] = 0 (108a) 

^ ^C^ll C02 ~ *^12 ^®1^®2 ^ *^13 ^®1^®2^®1 ^”*^*2 *^22 ^®1^®2^®1 

4r ^ ? 9 

- J33 C0-|S0^S02 - J-|3 50^002 - J23 S02(c^01 - s^0-|)]02 - (J33 C0-|S02 

*^13 ^®2 ~ *^23 ^®1^®2^®3 ” ^®2 ^2 ^®1^®2 ” ^3 C0^S02) 

+ [4 j*^ S02 C02 + 4 J*2 S0i(c^02 - 5^02) - 4 J22 s^0iS02C02 

- 4 J*3 c 0 ](c 202 - s^02) + 8 J23 S0TC01S02C02 - 4 J33 C^0^S02C02] 

- ^C'^22 *^33 2 *^23 ^®i^®i 3®2 “ ^*^12 ^®1 ^ *^13 ®®1^®1 

+ (J23 C0] + J33 S0])03 + K2 C01 + K3 S0] + fi C~'^12 ^®1^®2 

- J22 S01C01C02 - 0*3 S0^S02 + J23 " s^0t)c02 

+ J33 s0iC0]C02]> = 0 (108b) 
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■ '^*2 “l ^^*2 ^ ^^^^^22 ■ ■*■ ^^23 ^®1^^1®2 

+ 3^3 + J *3 c0^e2^ + K* e^ + K-j 00^03 

+ {[-2 J-|2 S02 ■*■ ('^1] “ '^22^ S0-]C02 ■*■ J23 c 0 iC 02 ] 9 i 

+ [(J 22 ■ ^* 1 ) C01S02 - 2J*2 s0iC0-|C02 + ^23 ^®1^®2 
+ J*3 (c^e^ - s^e^)c02]e2 + (J23 ^®2 " '^*3 501002)03 + s02 

- K-| s0iC02} "*■ ^ ^”'^*2 ^5^02 ” 3c^02) ~ J22)50"|^02^®2 

ic 0 0 0 ★ 

+ J-|2 s'^O'i (c^02 - 3 s 82) + 4J23 C0]S02C02 

- J^3 s9iC0t(c^02 - 3s^ 02)] - [-J*3 ©i + (J*3 C0i + 3^3 S0-|)02 

+ J33 03 + K3 + n (“^^*3 502 " '^23 501 C 02 + J33 C01C02)] ~ 0 ( 108 c) 

Considering Eq. ( 36 ) in conjunction with Eqs. ( 69 ), ( 99 ), ( 100 ), 
and (102), and letting i = 1 , we obtain the differential equation for v-j 

*^1 ^^^xl ^ ^1^ [(coiCa + a)2Sct) (ciJi Sa - u)2Ca) - 3fi^ (£g^-|Ca + £g2Soi) (Ag-jSa 
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o 2\ 2 
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• 9 

- £g 2^“)3 ''^1 (cotco-j + Sao)2)} “ '^'^1 ■*■ b^i)(uJiCo + <d 25 <*) 

+ {h^-^ + x^) 013 + v^] 

_ 2 2 

- ^ {(Jlcis« - i^cZ^ar + tj .3 + 2 - 3[(£giSa - £g 2 ^“) 
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(109a) 
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* 7 ^’z1 ''i'^i'^)* ;7 f-"zl 'i' '’ ■ I * Py1 ' ° • 


which is subject to the boundary conditions 


0 < x-| < Zi 


v^(0) =0 , vj(0) = 0 (109b) 
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Similarly, the differential equation for w.j is 
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S3 


+ (h^i + w^) [ 0)2 + 0)2 + 2n2 . 3 jj2 ^ 

- V-| (o)-|Co + o)2Sa)} " P] ■*■ ) (o)-|Sa - o)2Co) 

+ (hyi + V^)(o)^Ca + o) 2 Sa) + w^]_ ^ ^ 

+ £c 3 + 2 - 3 [(&g^sa - Ag 2 Ca)^ + ■*■ ^ 1 ^^ 

- + A^)} wj + f Elyi 
2 

+ “V ’ I Pzl ’ 0 » 0 < x^ < (llOa) 

9Xi ^ 

where w^ must satisfy the boundary conditions 

wi(0) = 0 , wj(0) = 0 (nOb) 

*”l '*’ [-o) 3 (o)iCa + o) 2 Sa) + 3 n 1^2 ''' 

■*■ (hyi + V.J ) [ 0 ) 3 ( 0 )-] So - o)2Co) - 3fi )ig 2 (^al^° “ ^a2^“^^ 

+ (h^T + w^) [u)i + o)| + 2^2 . 3^2 

• 3 

- V-] (o)-iCo + 0)2So) - — [(h^^ •♦• £]) (loiSo - 0 > 2 Ca) + (hyi + ) (o)-|Ca 


+ o)2So) + w^] - ■’■ *^c3 ■*■ 2 ■ 3[(«.g^sa 

- + ‘aal) C’xl + ‘1 )«i ) * I Ely, wjw;'^ 


9 Xi 


[EIy,W^' 


(1 - 




£i = 0 


EIy,«i' (1 



X, = I, = 0 




(llOc) 
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In the above 


co^ = - n(s62C03 + 591062593) - 06361 + O01S6262 

0)2 = ^{562593 - 591002063) + 563^ + 06106392 (Ilia) 

0)2 = fic0iC02 S0*]02 “*■ 63 

= C02C02 “ S0-|S02S02 

tg2 " -(062563 + 591S62O03) (111b) 

^a3 “ ^01502 

£^1 = -(502063 + 561062503) 

£o 2 = 502503 - 501002003 (lHo) 

to3 = 001002 


The equation5 of motion and boundary oondition5 a550oiated with the boom5 
2, 3, and 4 are obtained from Eq5. (109) and (110) by replaoing a by 
TT-cx, ir+oty and 2ir-a, re5peotively , and, of oour5e, ohanging the 5ub- 
5oript5 of vi , and Wi acoordingly. 

Following the 5ame prooedure a5 that U5ed to obtain Eq5. (109) 
and (110), the equation of motion and boundary oondition5 for Vg are 

P5 ^(^5 + ^5^ ^ -o)i(cj256 + U 3O6) + 3n^£ji(ta25® ^a3*^®^^ 

+ (hy5 + V5)[o)2 + U3 + 2n^ - 3 q^ (tj + ^a3^^ W5)[«1 (0)356 
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2 

- a) 2 C 6 ) - 3fi ( 1^3356 - (“ 3 CB + u) 2 S 6 )w 5 ) 

- P 5 [(\s + X5 )(o) 2C6 - 0)35$) - (h^g + W5 )(u 2SB + U3C6) + Vg] 

- <‘cl * <'c 2'8 - * 2 - 3[tJ + (tj 2 ce - «a35e)^]) 

<J»5 t(\5 ♦ «5)^ - ("x5 * ''5 * I E>z 5 ''5''5^} 

* 4 f-"z5 '5 (’ ■ f " Py5 - 0 • 0 < X5 < Z5 (112a) 

3X5 


VgCO) = 0 , v^(0) = 0 


I ”z5 vS' * 4 ^“^5 *5 • I “s'” 


^ ' ‘5 


= 0 


E«z5 ’B <' - 7 


^5 ” *”5 


= 0 


{112b) 


(112c) 


and those for Wg are 

P5 {(h ^5 + X 5 ) [((O 3 S 0 - a)2CS)(u2^^ “3^®) ~ ^^a2^^ 

O 

+ * V5)[a)i(o)2C6 - W3S6) - 3 n £gi (£92^6 - )^a 3 ^^^ 

< 2 2 2 2 2 
(i)] + (o) 2S3 + U3C3) + 2fi - 3fi [jig-j + (Jlg2^® 

- Vg (b)2S3 + W3C3) } 

- If P5 [-(h^5 + X5)a>^ + (hyg + V5 )(o 33C3 + a)2SS)] 

- jj- (*c1 * (‘‘c2‘^« - 2 - SC^ai * (»a2<ie ' *a3^8>^l> 

5 ' 

(1 P5 [(h,5 ^ »S)^ - (hxB * «b)^]> «B ^ | ^lyS "b^^) 


56 



ax 


2 t-^^y5'"5 ■ I Pz5 " ^ 


0 < Xc < «.t 


(113a) 


wJO) = 0 , w:(0) = 0 


(113b) 




Xg = ^5 


= 0 


Elyses 0 - f W5^)' 


(113c) 


^5 = ^5 


= 0 


The equations for boom 6 are obtained from Eqs. (112) and (113) by an 
appropriate change in subscript, and by replacing B by tt+b. 

b. Perturbation solution of the equilibrium problem . 

The first problem in attempting a solution of the equations of 
motion, Eqs. (108), (109), (112), and (113), is to identify the equilibrium 
configurations. To this end, we must let all the velocities and accelera- 
tions equal to zero in these equations. This leaves us three transcendental 
equations for the rotations ej (j = 1,2,3) and twelve nonlinear differential 
equations for the elastic displacements v^., w^. (i - 1,2 ,. .., 6). 

We shall consider the solution of the nonlinear equilibrium problem 
in the form 


Vio(Xi) - v.jgQ(x^. ) + v.Q^(x^.), w^q(x^-) = w^.qq(x^-) + w^g^(x^.) 


i = 1 » 2 , • • « 56 


(114) 


where the third subscripts on the right side of Eqs. (114) indicate the 
solution of the linearized problem if the subscript is zero and relatively 
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small perturbations if the subscript is one. It follows that the inertia 
matrix of the deformed body can be written as 

6 6 

[j( 0 )] = [otoiL + = E ]„[»,] + E [it.-fCOJlCHi] (115) 

U I 1 10' i=l ' 1 ' ' 

where is the inertia matrix as if the body was entirely rigid, in 

which 


'^0110 ' ^ 0 ’ '^0220 ^ ® 0 * *^0330 ^ ^0 
^^0120 ' '^0210 * '^0130 ^ *^0310 " '^0230 ^ '^0320 
are the moments of inertia of the rigid hub, and 


(116) 


'^illO 


^^1220 


'i330 


r^i ..2 .,2 


J 


2 ..2 


Pi(hy,.+h^^-)dx.+m.(hyi+hzi) 




0 • 
l 
0 
fi 
0 




' p j [( h^/x, IdXjtm, [ ( hj(, )^+1'yi 3 


= 1,2,... ,6 (117) 


'^1120 " ‘^iZlO 


rl\ 

0 

Ha 


"^ilSO 


= J 


1310 jQ 

rl 


Pi (h^i+x. )hy.dx.-m. (h^^.+^i )h^ • 

Pi(h^i+Xi)hzidx.-mi(h^i+£i)hzi 


J-ioon “ J.'< 


,-h . h_ .dx,--m. h...h 


1230 - ^1320 ■ 'Jq "iV'zl"^--'"l"yi zi 


are the moments of inertia of the appendages when in undeformed state, 
expressed In terms of local coordinates. Moreover, is the change 
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in the inertia matrix due to first order elastic displacements, which has 
the elements 


+ ni,(2hy,v,„o+Vj(,o+2h2,w,(,jtw^(,l,)^ _ 


'^i221 


‘^il31 


rH • 
fil * 

" '^1311 ' -Jd -"1(^i*>‘1)*'i00|x, 
“ ■^1321 ' -L'’1<\l*'i00*'’zi''l00*’l00“l00>'”‘l 


■”'l*''yi"iOO*''zl''iOO*''lOO"iOo) x,- = t, 

i = 1 , 2 ,... ,6 (118) 

To linearize the algebraic equations for the angles 9^ (j = 1,2,3) 
we would have to assume that the angles are small. This, however, is not 
always true for an arbitrary satellite, so that linearization cannot be 
justified. Fortunately, it is not difficult to solve the nonlinear 
algebraic equations for the angles ej (j = 1,2,3) by means of Newton- 
Raphson method for the moments of inertia given. As a first iteration. 
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we insert the moments of inertia of the satellite regarded as rigid into 
the three transcendental equations for 9. (j = 1,2,3), and obtain some 

J 

preliminary values for these angles. Hence, letting all terms in Eos. 

(108) involving time derivatives equal to zero, we obtain 

4J*2oC0iS02C02 + (J 220 " " 3s^02) + 43*3oS®l 

+ J23o(c^01 - s^9i)(3s^ 02 - c^e-j ) = 0 (119a) 

^ tItP 

4JiioS02C02 + 4Oi2qS0i(c ©2 - s^ 02 ) - 4O22qS^01 S 02 C 02 - 4J-| ^qC©] (c ©2 

- $^ 62 ) + 8J23gSe^ceiS92C02 - 4J22gC^9 -] $62062 = 0 (119b) 

-J*2o(®^®2 - 30 ^ 62 ) + 4(J*iq - J22o)s0iS 92C62 + J*2oS^9l (c^©2 - 3s^02) 

+ 4J22 qC 6 ^s 92C02 - J*2gS9] c9-| (0^62 - 3s^92) = 0 (119c) 

where [J*]. = [e]J [j(°ho [eJj in which [ejj is given by the first of 
Eqs. (104). Regarding the angles 6j (j = 1,2,3) as known constants, we 
can linearize Eqs. (109), (110), (112), and (113) with respect to v^., w.j 
and their derivatives, and solve for the perturbed elastic displacements. 
Hence, inserting Eqs. (114) into (109), we obtain the equations for v^j^^ 
in the form 

2 

Pi x-|) [((OiqCcx + o)2oSa)(«»)io^“ " ‘^20*"“^ " ^^alO^“ 

+ 2fi^ - 3n^^(£al0^“ ^a20^“^^ ^-aSO/O + (>i2i 
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- “20^^) • ^aSO^^alOS^ " ^320^“)]) - ^loo ^1 ^^l 

+ x-| ) 2^ ' *"020^“^ '*' ^c30 ^ ^ ~ 

^ ''a 30 ^> ''i'oo •^(^^ + M)^ - (^^ ^ 

+ m^(h^T + £l)} {(icioS“ - ^c20^“^^ ■" ^c20 ^ - 3[(^aioS“ 

- ^-a20^“)^ ■'■ ^a30^^ ‘ ''lOO (120a) 


where v-|qq is subject to the boundary conditions 

Vioo(O) = 0 , v'qq(O) = 0 (120b) 

2 

n^l {(*^xl '*’ ^1^ [(“10^“ ''’ “20^^^ ^‘^lO^^ ~ “20^“^ (^alO^“ 

2 

+ ^a20^“^^^al0^“ " ^a20'^®^^ ^'^yl ''l 00^1^ ^“10*^“ u>2QSa) 

+ o)§o + 2fi2 . 3o^ ((AalO^“ ^a20^“^^ ^a§o) 3 + (^i 

2 

+ Wiqq) [u)3Q(a)-]QSo - u)2oCo) - 3fi ia30(*'al0^“ " *'a20^“3] 


• _ 

- v-|QQ q2 + £^) ^(£^^qSo - £^20^®^^ ■'■ ^’cSO ■*■ ^ ■ ^I^^^'alO^® 

- ^a20Ca)^+ ^a3o) > + El v. 


Ml 

100 


f-T " 

''100 


Xi - 


Xi = £i 


= 0 


= 0 


(120c) 


The quantities ujo* HgjO ^cjO ^ “ 1»2»3) appearing in Eqs. (120) are to 
be calculated by using 0j (j = 1,2,3) as given by Eqs. (119). Note that 
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now primes designate total derivatives with respect to because v-|qq 
depends on x alone and not on t. Similarly, for Wiqq we have 

Pi {(hj^i + x-| )[- o)3o(wiqCoi + o)2oSa) + 3n^Aa3o(tgiQCa + 

+ (hyi + vioq) t‘^30^“10^°‘ ■ “20*-®^ ■ ^^^^aSO^^'alO^® " ^a20^“^^ 

+ (h^i + w-|qq) [(w^q + u)2o) + + Ag20^^^ 

- PI '^lOO (^1 + ((^cios® - ^c 2QC“)^ + ^c30 2 - 3[(£aio^® 

- ^ 320 ^“)^ + + Woo 4 PI - (^l + ^l)^^ 

+ m^(h^^ + Ai)} {(tj,^oSa " ^c20^®^^ * ^c30 * 2 - 3[(tgio^® 


- ^ 320 ^“)^ + ^a 30 ^^ - Woo = ° 

where w^qq is subject to the boundary conditions 

W]oo^O^ = 0 , w^'qq(O) = 0 

m, ((\, t -o)3o(uioCa + ti)2gSa) + 3fi'^tg3( tgigCo + ^a20^®^^ 

+ (\l ■'■ ^ 100 ^ r“ 30 ^“ 10 ^“ ■ ‘^ 20 ^®) " 3 fi‘^«.a 3 (*'al 05 ® " ^a 20 ^®)^ 

+ (h^i + w-|qq) [oj^q + 0)20 2^2 - 3^2 + tgio^^ 

■ Woo ^^xl ^ 1 ^ ( (^clO^® ■ ^c 20 <^®)^ ^c 30 2 - 3 [(£a-|o^® 

- ^a 20 ^®)^ ^ ^a§ 0 ])> Woo|x^ = M ^ 


(121a) 


(121b) 


(121c) 


Woo 


= 


= 0 
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The differential equations and boundary conditions fo,r and 
(i = 2,3,4) are obtained from Eqs. (120) and (121) by replacing the sub- 
scripts of and by the appropriate ones and the angle a by ir-a, 

TT+a, and 2TT-a, respectively. On the other hand, the differential 
equations and boundary conditions for v^qq and are obtained from 
Eqs. (112) and (113) in the form 

^^y5 ''500^ ^“20 “30 ^^a20 ^ ^aL^^ 

'^500^ t“10^“30^^ " “20^®^ ~ ^^^^al0^^a30^® ‘ 

“ V500 P5 (\s + ^5) ^^clO ^^c 2QC6 - <^c30^S)^ 

+ 2 - 3[£gfo + (£920^8 - £a3oSB)^]} + V500 {1. 

+ £ 5 )^ - (hx5 ^^clO (^c20^^ ■ ^c30^^^^ ^ 

O o 

- 3[£g^Q + (£320^8 - ^a30^®^ ''500 ^ ° (122a) 

where VgQQ satisfies the boundary conditions 

''500^0^ = 0 * ''^00^°^ = ° 

EIz5''5b0 =° • ^^z5''500 =0 (122c) 

as well as 

*^5 ^5^ ^(<^30^® ‘ o)2qC8)(u2oS 8 +. U3gce) - 3n^(£g3QS8 
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''SOO^ ^“ 10 ^^ 20 *^® ’ 

- 3ft^«,glo(«.a2oC6 - ■*■ ^^z5 ''■ '^500^ ^“lO (“20^® ^^30^®^^ 

+ 2 fi 2 - 3q2 [^Jq + (£^ 2 qS 6 + ^a 3 oC 6 )^]} - w^qq ogo^ (h^g 

+ X 5 ) a^io ^ (<^c 20 « - ^c30^6)^ + 2 - 3[?,Jo ^ (?.^ 2 qC£ - P. 33 oSe)^]) 

'*'500 '' 5 ^^ ■ (^x 5 ^ ^^clO ^ (^c 20 ^® 

- Mose)^ + 2 - 3[^afo ^ (^a 20 ce - ^a30SS)^3> - El^g wJJq = 0 (123a) 

" 500 ^ 0 ^-° ’ " 500 ^°^=° 


'^5 ^500 


El 


w: 


y5 500 


= 0 


(123c) 


"5 “ ^5 


The differential equations and boundary conditions for VgQQ and WgQQ are 
obtained by replacing In Eqs. ( 122 ) and ( 123 ) VgQQ and Wjqq by VgQQ and 
Wgoo and e by ir+ 6 , respectively. 

On the other hand, the boundary-value problem for the perturbation 
v^Q^ Is defined by the differential equation 

2 2 2 2 2 
Pl^ioi ^^“ 10 ^^“ “ 20 ^“^ “30 ■'■ n^alO^® ^a 20^“3 ''■ ^a 30 ^^ 

+ PlWioi [u 3 o(wtqSo - a> 2 oCa) - 3 fi A 33 Q(£aioSa ’ ^a 20 ^“^^ 

+ v^Q^ (h^^ + x^) ^(^cl0S“ “ ^020*^“) + ^c30 ^ 
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- 3C(^,o=» ■ 'a2o"“>^ * \lo^) * ''i'o0'''i'00 

* ''i'oi - (^'xl * '"l(^l 

* m|> ( (‘clO^'‘ ■ ’■C20'=“> * ‘c30 2 ■ 3[(«a,o^'’ ' 

* ‘alol) * EIjl Ct (''i'00>^ * ''i0O’'l'ooJ> 


5''i'oo’'ioo]* 


’•a20C”) 


2 


* Ela, ''I'oi ''ioo'-ioo) ■ Elzl ''iol 

* ''ioo''i'oo'''i'oo * I ''ioo 


f’ - 1 (''ioo)^] ' -E'zitf(vi'oo)^ 

(124a) 


and the boundary conditions 


vioi(O) = 0 > ''101^°^ " ° 

"'l"l01 C“30<"10®“ ■ “20'=“) - ^“^S30<‘‘al0*“ ' *a2o'“)) 

♦ vjoi {m,a2 (h^, t «,) ((t^iQSa - t^jgca) ♦ tJo 2 
- 3[(*al0*“ - "320“)^ * *a30)> ‘ | ^Iz, [2v;„(,v'foo a 

■ ''loi ^®''ioo''ioo^ 

^ ''i'bi “ I * I ^^zi ''Too ) 




- 1 (''loo^^^ ''i'oi 

®''ioo''ioo''ioi^ ' ■ l^^zi ^''1oo^^''ioo 


V at (124c) 
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and that for the perturbation by 

PlVioi ^a30^*'al0^“ “ ^a20^“^^ 

+ Piw^Q^ [uifo + <^20 + 2n^ - 3^2 (ta?0 + ^alo^^ 

+ wjoi {-Pi«^ {\i + Xl) ((^cl0S« - tc20<^ot) + Ac3o + 2 

- 3[(iial0S“ - ^a20^®^^ ^a30^) ^^yl ’^lOO'^lOO 

+ 5 y^'iMooh + w^'oi pT [(h^l + Ai) 2 . (h^i + xi)2] 

+ m■^{h^■^ + {(^clO^® ■ ^c20^“) ^c30 * ^ 

- J^a20^“)^ ^a30^) ^^yl Hoo)^ ■*■ '^lOO'^lOO^^ 

+ ElylW^'Ql (10 WjQQW^Qg) 

“ '^101 ■ f “ ‘ ^^yl *-f '^loo'^ioo'^l'oo 

w^Q^(0) = 0 . ° 

2 

"’I'^IOl ^“30^“10^“ " “20^“^ “ ^a30^*'al0^“ ~ *’a20^“^^ 

+ [oj^Q + 0)20 * ^^alO * ^a20^^ 

+ wjg^ (h^^ + { (^clO^® ■ ^c20^“^ ■*■ *^c30 '*’ ^ " 

- ^a20^“)^ ^ ^a303> ’ f ^lyi C2wi'ooWl”oo ^ 
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- Ely! w\'oi (5w-jQQwyQQ) 

* Elyi [1 - I (wjoo)^] " f Ely, [(wioo>^ « 

+ 2wjQ0 (wj'Qg)^] 


Ml 

100 


X5 = ^5 


Ely, to - I (Wio)^] 

''100”! 00”! 01 ■' " ■ 2 “yl'”100^ "100 


] = - ^ EI„i (w-J«^) wVf 


X5 = ^5 


(125c) 


with companion equations for v^*qi and w^qi (i = 2,3,4). In a like manner 

P5'^501 I-“20 “30 ^ ~ ^^a20 * ^a30^^ 

■*■ P5'^501 [^“10^“30*® " “20^®^ “ ^^^^al0^^a30^® " ^320^^®^^ 

^501 {-PsO^ (hj(5 + Xg) + («-c20^® ■ ^c30*®^^ ^ 

- 3[^al0 + (WS - ^a30S6)^]) + E^zS + 5 v^5qV^qO^> 


^ ''501 


(} P5 L(hx5 + ^5)^ - (hx5 + xs)^]^ ^ £^.^0 + (^c20^^ 

- ^c30S®)^ 2 - 3[tjQ + (Aa20^^ ' ^a30^®^^^^ ■*■ ^''500^^ 

" 10 ''500''5W> " E1z 5 ''501 ^1° ''Wsoo) 


- 'Iz5 ''501 n - I (v^oo)'3 = - [f (v^oo)' 
10 ''500'' 500'' 500 ^ f ''500 ^''500^^^ 


(126a) 


''501^°^ ° ’ ''501-^°^ ^ 


(126b) 


67 


f ''500''500 ^''500^^^ ''501 " ^ ''500''500''501 


^n-f(v^oo)2] v^5^) 

" I •^(''500^^ ''500 ^ ^ ''500 ^''500^^^ i 

i 

■ ^^z 5 ■ f ^'' 500 ^^^ ''501 I ^5 " ^5 

“ ^ '' 500 '' 500''501 ^ " ■ I ^^z 5 ^'' 500 ^^ ''bOO_ 


and 

OgVgo') C‘^o^‘^0^^ ■ “30^®) ■ 3f^^alo( ^20*^^ ■ ^aSO^s)] 


^ pb'^ 501 ^“^lO ('^20^® ^^alO "*" ^^a20^® 

^aSO'^s)^]) + {- P5fi^ (hj^5 + Xg) + (^^20*^^ 


- ic3oS6)^ + 2 - 3[£a?o (^a20C6 " ^^a30S6)^]) 

+ Elyg [10 W500W500 + 5 W5ooW5oo]> 


+ WgQ^ {^^Q-Ps C0ix5 ^5)^ - (hx5 ^5)^]) ( ^clO (<^c20^e 
- ^30^®)^ ^ “ ^l^^alO ^ ■ ^a30^®)^0 ^^5 ^'"bOO 

+ 10 WgQQWg'oo]} + Elyg w^"oi '^500''^500^ 


+ 


^^yB '^501 ’ f ^'^BOO^^^ ~ “ ^^y5 

I II III 5 **** /* \?T 

’^BOO'^BOO'^BOO ^'^500 ^'^500^ ^ 


[| (w 



3 


Wbo-|( 0) - 0 , w^Ql(0) = 0 


(126c) 


2 

(127a) 

(127b) 
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f ” 500'^500 ^*^ 500 ^^^ '^501 “ ^ '^ 500 '^ 500'^501 


+ [1 - I (w5oq)^3 

^ I ^^y5 '^500 ^ '^500 ^'^500^^^ 

- Elyg {[(1 - I (w5oq)^] w^'oi 


i at Xg = £5 ( 127 c) 


■ ^ '^500'^500'^501 > ^ ■ 2 ^^y5 Koo^^ '^500 


It should be pointed out that this particular perturbation scheme 
enables us to solve first Eqs. (119) for the first approximation rotations 
9jOo (j = 1*2,3) independently of the elastic displacements. The rotations 
are then introduced into Eqs. (120) - (123), yielding the first approxi- 
mation for the elastic displacements v^-qq and w.|qq (i = 1,2,..., 6) 
independently of one another. Inserting the first approximation v^qq 
and w^-QQ into Eqs. (124) - (127), we can obtain the corrections v^-q-| and 
Wioi to the elastic displacements. The sums of these solutions yield v^-q 
and w^-Q (i = 1,2,. ..,6) according to Eqs. (114). Then, inserting v^-q and 
'"^iO * 1 *2,. . . ,6) back into Eqs. (119), we obtain the angles ejQ 

(3 =1,2,3). In the vast majority of cases, this approximation is suf- 
ficient. If not, having the new angles, we can iterate once more to 
improve the elastic displacements v^-q and w^-q, as well as the angles Gjq. 

c. Liapunov stability analysis and the eigenvalue problem . 

The values Sjq, v^.q, and w,-q obtained above, together with sets 
of admissible functions 'J'j(x^-) and i|»j(x^-)» are subsequently introduced 
into Eqs. (49), (51), and (53), to obtain the coefficients m^j^, fj^* and 
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kji^. The coefficients nij|^ and yield the symmetric matrices [m] and 
[k], whereas using Eq. (57) the coefficients fj|^ yield the skew symmetric 
matrix [g]. 

From Sec. 6, if [k] represents a positive definite matrix, then 
the nontrivial equilibrium is asympototically stable. On the other hand, 
to obtain the natural frequencies, we must solve the eigenvalue problem 
in the form (66). However, before the nontrivial equilibrium can be 
determined, the system stability tested, and the natural frequencies cal- 
culated, it is desirable to use specific values for the system parameters. 
This is done in the next section. 


d. The shortening of the projections effect 

As indicated in Sec. 2, the booms are assumed to be inextensional , 
so that there is no longitudinal vibration. However, because of the 
transverse displacements, there is a shortening of the projection on the 
nominal axis of any element of length of the boom. In fact, from 
Eq. (14), the change in length of projection of any element of length dx^ 
is 


= _ 1 


fav/ 

2 

1 

3W^. 

2- 


+ 


2 




9X,- 

1 N 

1 


dx. 


i = 1 ,2,. . . ,6 


(128) 


We shall treat this shortening as a perturbation of the spatial 
coordinate x^- , so that we can write 

^i " ^iO ^il ’ • i = 1,2,...,6 (129) 

where x^g are the original spatial coordinates and x^.^ are the 
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perturbations. From Eqs. (128), however, we conclude that the shortening 
is a second-order effect. Hence, it will not affect Eqs. (120) - (123) 
except that x^. are to be regarded in these equations as x^.q (i = 1,2,..., 6). 
This enables us to solve for v^-qq and w^^q and write the shortening of, 
the projections in the form 



r* 

f «> 

2 

f y 

2-, 



S^iOO 

4. 

^'^iOO 


- 


X 




i = 1,2,... ,6 (130) 

where is a dumny variable. On the other hand, the perturbation equa- 
tions, Eqs. (124) - (127), must be modified to account for the shortening 
effect. For example, the boundary-value problem for v^-qi becomes 

PlVioi {((*>ioCa + o)2oSa)^ + «30 ^a30^^ 

0\ 

+ Pl^ioi f“3oHo^“ ■ “20^“^ ■ 3^^*a30 ^^-alO^® “ ^a20^°^^ 

+ vjQ, {- (h^^ + x^q) ((^c10^“ ■ ^c 2QC°^ ^c§0 + 2 

■ ■ ^a20^“^^ ^aSo)'*’ ^^zl ''lOO''lbo ^ ''iWoO^^ 

''101 ^ ^ "’l^^l 

^10^} (^\l0^“ ■ ^c20^“^^ ^c§0 ^ 2 - 3[(ilaioSa - ^a20^“)^ ^aSO^,) 

^ ^^zi Ct (''i’oo)^ ^0 ''ioo''ibo^> " ^^zi ^vroi no vioovfoo) 

- ''iol n - I (v^oq)^]} = -EI,i [f (vYoo)^ - 10 v^00''V00Vl00 

I ''100 ^''ioO^^^ ^ *'1 nuiioCa + a)2oSa)(o)^QSa - co^QCa) 
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2 




^10 P 


0 


3v 


100 


Hi 


uV 


qT' 

0 

0 

2' 

9^1 J 



dC] , 0 £ x-|Q <_ £ 


10 


(131) 


subject to boundary conditions (125b) and (125c), where in the latter 
£] must be replaced by the shortened length £^q. 

12. The RAE/B Satellite. Numerical Results . 

The general formulation of Sec. 11 has been used to obtain the 
nontrivial equilibrium configuration of the RAE/B satellite, to test the 
stability of equilibrium, and to calculate the natural frequencies of 
oscillation about the nontrivial equilibrium. The system parameters are 
as follows: 

Ag = 87.74 slug ft^. Bg = 83.74 slug ft^, Cq = 18 slug ft^ 
p] = P 2 = P 3 = P 4 = 4.348 10“^ slug ft"^ P5 “ P 6 ' ^*596 

10 "^ slug ft“^ 

m-| = m2 = m3 = m^ = 2.40 ^ 10"^ slug, mg = mg = 0 
£] = ^2 “ ^3 “ ^4 ~ 600 ft, £5 = £g = 315 ft 


• — 1 
UJ 

= E'zl 

- EIy2 

= 

• • “ 

EI,4 = 

15.278 lb ft^, a = 

30“ 


= EIz5 

= Elye 

= EI26 = 

13.889 lb ft^, e = 25“ 


^xl * 

^x4 

0.973 

ft, 

^2 = 

•^x3 “ 

0.878 ft, hj(5 » hjjg 

' 0 

hyl = 

"V ' 

to 

0 

0 

ft, 

•^y2 

' "^y3 

* -0.760 ft, hyg ■ 

^y6 

^1 * 

h^2 - 

^z3 “ ' 

^4 



- 0 


q ■ 

4.653 

X 10"^ 

rad 

sec" 

1 




We shall present the results of the analyses in the order listed above. 
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a. Nontrivial equilibrium 


Inserting the above data into Eqs. (116) and (117), as indicated 
in Eq. (115), and solving Eqs. (119), we obtain 

0100 " 0.13537 rad = 7.756 deg. 

0200 = -5.63789 x lO"^ rad = -3.2302 x lO"^ deg. 

0300 = 1.37374 x 10“® rad = 7.87096 x lO"® deg. 

We note that 0.|oo is caused largely by the damper booms. The fact that 
the rods are not attached at the satellite mass center turns out to have 
an insignificant effect on 0joo (j = 1.2,3). 

To evaluate the elastic displacements v.|oo(Xt) and w.jQQ(x.j) 

(i = 1,2,..., 6). We assume the solution of Eqs. (120) - (123) in the form 

''iOO^^i^ “ ®ri0^r^^i^ 

i = 1,2,. ..,6 (132) 

P 

r= I 

where 

<l>y.(Xj) = A^[(cos + cosh 6^Jl.j)(sin B^x.,- - sinh B^Xj) 

- (sin + sinh BpAi)(cos B^x^- - cosh B^Xj)] (133) 

are eigenfunctions corresponding to a bar in bending with the end x.j = 0 
fixed and having a mass m^ attached at the end x^ = The eigenvalues 
are solutions of the characteristic equation 
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(134) 


mi 

(1 + cos cosh BpA-j) = By.iii* — — (sin cosh B^A-j 

^i*'i 

- cos sinh B^a-,*) 

Moreover, the amplitudes A^. are such that the eigenfunctions '{'^(x^*) are 
orthonormal , i.e., they satisfy the relation 


f^i 

j P'1 ^)p(x^ ) <}>g (xj )dx.| "t" m.j (j)y,( ) <{ig ( ) “ *^irs (135) 


where is the Kronecker delta. Limiting the series in (132) to two 
terms, p = 2, the first two roots of Eq. (134) and the amplitudes A^. 
corresponding to i - 1,2, 3, 4 are 

B^t^- = 1.85813 A^ = 0.47696 

B^Z. = 4.65310 Ag = 0.03789 

In addition, the coefficients a^^.Q, (i = 1,2, 3, 4) are 

Table I. 


i 

®li 0 

® 2 i 0 

*^110 

*^ 2 i 0 

1 

- 0.13656 X 10 ^ 

- 0.98055 X lO'"' 

0.55184 

0.46385 X 10 '^ 

2 

0.13652 X 10 ^ 

0.97981 X 10 “^ 

0.55188 

0.46385 X 10 “^ 

3 

- 0.13652 X 10 ^ 

- 0.97980 X 10 '^ 

- 0.55187 

- 0 . 4638 A X 10 '^ 

4 

0.13656 X 10 ^ 

0.98054 X 10 ’^ 

- 0.55184 

- 0.46385 X 10 "^ 


slug' 

slug' 


■ 1/2 

• 1/2 


( 136 ) 


The first two roots of Eq. (134) with = 0, and the amplitudes A-] and 
A 2 corresponding to i =5,6 are 
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= 1.87511 = 0.63510 slug“^/^ 

(137) 

g2ili = 4.69414 A2 = 0.04899 slug'^^^ 

whereas the coefficients a^^.Q and are 


Table II 


i 

®li0 

^2i0 

^liO 

^2i0 

5 

-0.93855 X 10'2 

-0.12900 X 10"^ 

0.17756 

0.70688 X 10'^ 

6 

-0.93840 X 10'^ 

-0.12900 X 10"^ 

0.17756 

0.70688 X 10"^ 


It will prove of interest to list the elastic displacements of the end 
points, as calculated by means of the linearized equations. These dis- 
placements are 


Vioo(i^l) = -52.205 ft. 

V2Qo(t2) = 52.192 ft. 

VSOO^^S) = -52-191 ft, 

V4oo(t4) “ 52.204 ft, 

''500^^5^ = - 4.8655 x lO"^ ft, 
V6 oo(J^ 6) = -4.8648 x 10"^ ft. 


Wioo(^l) — 2.1071 ft, 
'*^200^^2^ “ 2.1072 ft, 
Wsoo(t3) - “2.1072 ft, 
W4 oo(M) = -2.1071 ft, 
W5 oo(*- 5) = 0.92961 ft, 
'^600^'^e) = 0.92960 ft 


The above values of v^.qq(x^.) and w^.qq(x^-) (i = 1,2,... ,6) enable 
us to solve Eqs. (119) for the angles 6jg (j = 1,2,3) and Eq. (131) and 
the companion ones for the perturbations v^.q-|(x^), w,*gi(x^-) (i = 1,2,..., 6). 
The resulting angles are 
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P]Q = 0.1Q695 rad = 11 .?846 deg. 

= -6.54250 X 10*^ rad = -3.74858 x 10"^ deg. 
03 Q = 1.37608 X 10"^ rad = 7.88438 x 10"^ deg. 


Instead of listing the perturbations v.jqi and w.jqi, we shall list the 
complete solutions v.jq and w.jg in the form of the series 


r=1 


i = 1,2 6 


(138) 


Wio(Xi) = 


r=1 


where are still given by Eqs. (133), in which the eigenvalues 

and amplitudes (r = 1,2) are given by (136) and (137). The final 
results are tabulated as follows 


Table III 


i 

^li 



®2i 


hi 

hi 


1 

-0.13803 

X 

10^ 

-0.86379 X 

10“'' 

0.54865 

0.46378 X 

10-2 

2 

0.13800 

X 

10^ 

0.86313 X 

10'^ 

0.54869 

0.46378 X 

10"2 

3 

-0.13800 

X 

10^ 

-0.86313 X 

10-1 

-0.54869 

-0.46378 X 

10"^ 

4 

0.13803 

X 

10^ 

0.86379 X 

10-1 

-0.54865 

-0.46378 X 

10-2 

5 

-0.93855 

X 

10'2 

-0.12900 X 

10-^ 

0.17756 

0.70670 X 

10-^ 

6 

-0.93840 

X 

10"2 

-0.12900 X 

10-^ 

0.17756 

0.70670 X 

10-^ 
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Moreover, the final end di ^ pi acementn are 


= -52.816 ft, 

V2o(f'2) = 52.803 ft, 

V3o(t3) = -52.803 ft, 

V4q(*^ 4) ” 52.816 ft, 

V 5 o(f. 5 ) = - 4.8655 x 10"^ ft, 
VeoC^e) = - 4.8648 X 10"^ ft, 


w.,^ = 2.0948 ft, 

= 2.0950 ft, 
W3 o(^ 3) = -2.0950 ft, 
w^q(;. 4) = -2.0948 ft, 
Wso(t5) = 0.92962 ft, 
WgQ(£g) = 0.92962 ft. 


and we note that the nonlinear effect is virtually zero for booms 5 and 6. 
The nontrivial equilibrium is depicted in Fig. 7, where only the radial 
booms are shown because the displacements of the damper booms are in- 
significant. 


b. Liapunov stability analysis 

A stability analysis using k, as given by Eq. (60), as a testing 
function has been carried out. Essentially, the analysis reduced to 
testing the matrix [k] for positive definiteness, where the elements of 
[k] are given by Eqs. (53). The numerical values of the elements for 
the particular configuration at hand are listed in the next subsection. 
The matrix was found to be positive definite, so that the equilibrium is 
asymptotically stable . 

c. Eigenvalue problem 

Using Eqs. (49), (51), (53), and (57), in conjunction with the 
above data, we obtain the elements 
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3.42090 X 10 0 -6.31727 x 10"^ 4.78807 x 10"^ -7.67484 



-7.67484 X 10 -1.57707 



-4.78845 X 10 ' 3.44664 x 10 1.77064 
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-6.74965 X 1.14635 x 10 ^ -1.24054 



-6.74931 X lO'-^ 1.14630 X 10"^ -1.24288 



-3.36689 x 10'" 2.64824 x 10'-’ 1.40028 
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-7.49380 




-5.43408 X 10-6 1.72556 x 10-5 6.59170 
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-3.60230 



- 8.10040 X 10 3.88945 



Note that the elastic displacements were represented by one mode each. 

Using the formulation of Section 7, we obtain *the following natural 
frequencies 

Table VII 


+ 1.514583 X 10“^ 

+ 9.554415 X 10’^ 

+ 9.900792 X 10'^ 

, + 6.143208 X 10‘^ 

(ji. - rad sec" , 

^ + 4.319115 X 10"'^ 

+ 2.020401 X 10"^ 

+ 2.020401 X 10"^ 

+ 1.975175 X 10"^ 


+1.955743 X 10"^ 
+1.871461 X 10"^ 
+ 1.856504 X 10"^ 
+1.856504 X 10"^ 
+ 8.708815 X 10"^ 
+ 6.155318 X 10"^ 
+ 3.439481 X 10"^ 


The natural modes have also been obtained and will be discussed in a 
future paper. 

d. Parametric study 

The stability analysis was carried one step farther by varying the 
angle a. It was found that the system was asymptotically stable for 
a = 50°, but became unstable for a = 51°. The results can be easily 
explained by the fact that in the absence of damper booms and for com- 
pletely rigid radial booms the system becomes unstable around a = 45°. 

The gravitational and centrifugal effects tend to deform the flexible 
booms in a manner that the moments of inertia about the local vertical and 
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about an axis tangent to the orbit are the same for an angle a such that 
50® < 0 < 51®. It should be mentioned that instability in both cases can 
be traced to angle 03 , which tends to become large when the moment of 
inertia about the local vertical becomes larger than that about the axis 
tangent to the orbit, as at this point the "least moment of inertia" 
criterion is violated. 

The same parametric study was undertaken with respect to the 
natural frequencies. In terms of natural frequencies, instability occurs 
when at least one natural frequency (we recall that in our case the natural 
frequencies occur in pairs) reduces to zero. Here again the system be- 
comes unstable for 50® < a < 51®, thus corroborating the results obtained 
by the Liapunov stability analysis. 

13. Surrmary and Conclusions 

Two new theories for studying the motion characteristics of a 
rotating system with flexible parts about undeformed equilibrium have been 
developed. The first is qualitative and the second quantitative. 
Specifically, the first represents a stability theory and the second a 
method for obtaining the system natural frequencies. 

The stability theory is based on the Liapunov direct method and 
makes use of modal analysis to represent elastic displacements. The 
novelty of the formulation lies in the fact that for the first time a 
nontrivial equilibrium is considered in conjunction with the Liapunov 
direct method for a stability analysis of spinning flexible bodies capable 
of large deformations. 

The stability analysis can be divided into two major parts: the 

evaluation of the nontrivial equilibrium and the stability analysis itself. 
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When the body is capable of large deformations, nonlinear algebraic 
and differential equations must be solved for the rotational and elastic 
displacements, respecti vely, where these displacements define the 
equilibrium configurations of the system. Because the problem is one of 
stability about nontrivial equilibrium, it is necessary to expand the 
Liapunov function about that equilibrium. Assuming small displacements 
from equilibrium, the problem reduces to the evaluation of a Hessian 
matrix at the nontrivial equilibrium and testing the matrix for sign 
definiteness by means of Sylvester's criterion. It should be pointed out 
that the size of the Hessian matrix depends on the number of eigenfunctions 
used to represent the elastic displacements. 

The method for obtaining the natural frequencies of the system 
makes use of the variational equations about the nontrivial equilibrium. 
Then the set of second-order differential equations is converted into a 
set of twice the number of first-order differential equations. The 
associated eigenvalue problem yields the system natural frequencies. 

The two methods are quite general in scope, and can be used for 
testing stability and calculating the natural frequencies of a large 
variety of hybrid systems. As an application, the theory has been used to 
test the stability of the RAE/B satellite. First, the nonlinear equations 
have been solved for the nontrivial equilibrium configuration, and then 
this configuration has been used to evaluate the associated Hessian 
matrix. The satellite was found to be stable. Then one of the systems 
parameters has been varied to predict at which point the equilibrium be- 
comes unstable. The results are in line with the expectations. In 
addition, the system natural frequencies for oscillation about the 
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deformed equilibrium were calculated. The parametric study used in con- 
junction with the Liapunov stability analysis was used to examine how the 
frequencies are affected. The study resulted in the same instability 
statement. 
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FIGURE 5 - DEFORMATION OF ELASTIC BOOM 
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FIGURE 6 -RADIO ASTRONOMY EXPLORER - LUNAR (RAE/B) SATELLITE 


BOOM I 



FIGURE 7, NONTRIVIAL EQUILIBRIUM CONFIGURATION 
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